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ABSTRACT 

The Standard Model of particle physics makes it possible to simulate 
complete events for physics signatures and their backgrounds in high en- 
ergy collisions. Knowledge of how the produced particles interact with the 
materials in a detector makes it possible to simulate the response of any 
particular detector design to these events and so determine whether the 
detector could observe the signal. The combination of these techniques has 
played an important role in the design of new detectors, particularly those 
for hadron supercolliders where the high rates and small signal cross sec- 
tions make the experiments very difficult. The technique is reviewed here 
and illustrated using the simulations of the GEM detector proposed for the 
Superconducting Super Collider. Although the simulations and results de- 
scribed here are somewhat detector-specific, we believe that they can serve 
as a useful model for this component of detector design for future hadron 
supercolliders. 
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1. INTRODUCTION 

The "Standard Model" of particle physics describes strong, electromagnetic, and weak 
interactions in terms of an SU{3) (8) SU{2) U{1) gauge theory. The theoretical elements 
of the model have been in place for more than 20 years. ^'^ Over that period a series 
of experimental results — including the discoveries of charm,^ of the r-lepton and the b 
quark, and of the predicted VF^ and Z'^ bosons^ — have provided increasing evidence 
of the correctness of the Standard Model at present energy scales. More recently, the 
Standard Model has withstood extremely stringent quantitative experimental tests. ^ Thus 
it is now accepted that it provides a good description of particle physics down to the 
shortest distances yet probed, about 10~-'^^cm. 

In the Standard Model the basic constituents of matter are spin-| quarks and leptons. 
There are six types or "flavors" of each. Under the SU{2) ® U{1) group of electroweak 
interactions, the left-handed fields are in three quark and lepton doublets; the right-handed 
ones are in singlets: 

Lepton.: h\ i""'] f"'") (1.1a) 
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Under the SU (3) group of strong interactions, which for historical reasons is known as 
color, the quarks are triplets and the leptons are singlets. The quarks and leptons can also 
be organized into three generations identical except for mass. There is mixing between the 
generations of quarks; the primed lower elements in Eq. (1.16) denote linear combinations of 
the mass eigenstates. All the fermions definitely have been found except for the top quark 
and the tau-neutrino. Direct evidence for the top quark with mass mt = 174 ± 10 GeV 
has recently been reported.^ There is strong indirect evidence for the tau-neutrino.^'^ 
If the number of quark-lepton generations is equal to the number A^i^ of light neutrinos, 
then there are no more than these three. The evidence for this comes from precision 
measurements of the at LEP, which imply Ni, = 2.99 ± 0.04 in the Standard Model. ^ 

The interactions of the quarks and leptons are invariant under SU{3) (8) SU{2) (8) U{1) 
locally at each point in space-time. The requirement of local invariance implies that there 
is one spin-one gauge boson for each generator of the symmetry group, and it restricts their 
couplings so that the theory is renormalizable. This means that cross sections do not grow 
more rapidly than allowed by unitarity and that the form of the interactions is insensitive 
to possible physics at much higher mass scales. These gauge bosons are responsible for all 
known interactions except gravity, for which there is no fully satisfactory quantum theory. 
The eight generators of the SU{3) color group correspond to the eight massless gluons 
of Quantum Chromodynamics (QCD) that are responsible for strong interactions. Gauge 
invariance requires that these interact. These self-interactions produce a potential energy 
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which grows hnearly with distance between isolated quarks or gluons, which are therefore 
permanently confined into the observed hadrons such as the proton. At short distances — 
large energies and momentum transfers — the same interactions make the QCD coupling 
weak, so that perturbation theory can be used to describe the production of heavy particles 
or large-pT interactions in hadronic collisions. Because a gauge theory restricts the form of 
the couplings and because calculations can be done in perturbation theory, the Standard 
Model can be used to predict cross sections even for new physics (provided, of course, that 
a model exists for the couplings of the new physics to standard model degrees of freedom). 

The remaining four generators of SU{2) ® U{1) correspond to the massless photon 
and the massive and bosons responsible for electroweak interactions. If a gauge 
symmetry is unbroken, then the gauge bosons and fermions associated with it must be 
massless. It is possible, however, to introduce new interactions that break the symmetry 
"spontaneously" while preserving the good high-energy behavior of the gauge theory. This 
is known as the "Higgs mechanism" In the Standard Model, both the and Z'^ bosons 
and the fermions get their masses from an elementary scalar field, which is a doublet under 
SU (2) and which acquires a vacuum expectation value, thus breaking the symmetry. 

While there is strong experimental evidence supporting the gauge-theoretic part of 
the Standard Model, there is as yet no evidence for or against the Higgs mechanism for 
electroweak symmetry breaking. Nor is there data to indicate the mechanism by which 
finite and unequal fermion masses are generated (flavor symmetry breaking). Understand- 
ing how these masses are produced is the central problem in particle physics. Several 
scenarios have been proposed for electroweak and flavor symmetry breaking: 

• Standard Higgs models, containing one or more elementary Higgs boson multiplets. 
These are generally complex weak doublets. The minimal model has one doublet, with 
a single neutral boson H*^. 

• Supersymmetry. In the minimal supersymmetric standard model there are two Higgs 
doublets, and every known particle has a superpartner. 

• Models of dynamical electroweak and flavor symmetry breaking. The most studied 
proposal is technicolor-plus-extended-technicolor, with one doublet or one family of 
technifermions . 

• Composite models, in which quarks and leptons are built of more fundamental con- 
stituents. 

None of these proposals is fully satisfactory. In elementary Higgs boson mod- 
els, whether supersymmetric or not, there is no plausible natural explanation of why 
electroweak symmetry breaking occurs and why it has the observed scale. In non- 
supersymmetric models, the Higgs boson's mass, Mjj, and its vacuum expectation value, 
V, are unstable against radiative corrections. There is no natural reason why these two 
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parameters should be very much less than the energy scale at which the essential physics 
changes, e.g., a unification scale or the Planck scale. This instability does not arise in 
supersymmetric models, which is the motivation for supersymmetry at the electroweak 
scale. ^^'^^ Furthermore, elementary Higgs boson models are known to be "trivial", i.e., 
they cannot make sense as interacting field theories with the cutoff taken to infinity. 
This means that they must be regarded as effective theories, meaningful only below some 
cutoff A. Obviously, A must be somewhat greater than Mfj for the effective theory to 
make sense. For a modest separation of these energies, Mjj < few x A, both perturbative 
and lattice gauge calculations give Mjj < 650 GeV in the minimal one-doublet model. 
Finally, elementary Higgs models provide no clue to the meaning of flavor symmetry and 
the origin of its breaking. The flavor-symmetry breaking Yukawa couplings of the Higgs 
boson to fermions are arbitrary free parameters. 

Despite these apparent problems, the standard Higgs boson, H^, charged Higgses, H^, 
and the supersymmetric partners of all the known particles may exist and must be sought. 
However, if something like the standard H'^ is found and is heavier than about 700 GeV, 
experiments must have the capacity to discover the additional, unspecified new physics that 
surely exists in the same energy region, of order 1 TeV. One of the basic considerations 
driving the high collision energy of the SSC was the possibility that such physics exists. 

Dynamical theories of electroweak and fiavor symmetries — technicolor and extended 
technicolor — address these shortcomings of the elementary Higgs boson models. How- 
ever, they do so at the price of introducing flavor-changing neutral currents that are too 
large, and pseudo-Goldstone bosons (technipions) that are too light. Further, it is diffi- 
cult to build realistic models with QCD-like dynamics that are consistent with the precision 
tests of the electroweak interactions. These difficulties have been mitigated, but only by 
invoking a new form of strong technicolor dynamics. ^^'^"^ Realistic models of composite 
quarks and leptons are similarly difficult to construct. While no compelling models of 
dynamical electroweak and fiavor symmetry breaking exist, there are nevertheless model- 
independent phenomenological programs for searching for them.^^ 

While the origin of electroweak symmetry breaking is unknown, the mass scale for 
the new physics responsible for it is set by the vacuum expectation value of the Standard 
Model Higgs boson: 

1 _i 

v = 2-iGp^ = 246 GeV. 

The cross section for any new physics at this scale is set by dimensional analysis and 
coupling strengths to be cr ~ 1 nb - 1 fb. So long as the nature of this new physics 
remains unknown, the most promising way by far to explore this mass and cross section 
scale is with a pp (or pp) supercollider. That is because such a machine functions as a 
collider of broad-band beams of quarks and gluons (thus allowing coupling to both color 
and electroweak quantum numbers) and because these beams have the highest energy and 
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luminosity technically feasible today. Clearly, it would have been easier to carry out the 
high-mass searches at higher center of mass energy and a luminosity of £ ~ 10"^"^ cm^^ s~^ 
rather than 10"^"^ cm~^ s~^. Nevertheless, for the foreseeable future, the capabilities of the 
CERN Large Hadron CoUider (LHC), = 14TeV and £ = lO^^cm^^g-i^ j^^]^^ ^j^g 
best instrument for discovering the origin of electroweak symmetry breaking. 

Hadron supercolliders such as the SSC or LHC present a difficult experimental problem: 
the interesting cross sections are very small compared to the total cross section of order 
100 mb and small even compared to the dijet cross section at a comparable mass. Thus, 
careful study is needed to determine whether a given experimental design can detect the 
signatures of models for new physics. Because these signatures always manifest themselves 
as quarks and gluons (jets), charged leptons and neutrinos (missing transverse energy, 
^r), the Standard Model allows simulation of events both for the new physics and for its 
backgrounds in a reliable and predictive way. Given a detector design and known facts 
about particle-matter interactions, the response of the detector to these events can also 
be simulated. Finally, the simulated data can be reconstructed and analyzed just as for a 
real experiment to determine whether one could observe a signal. 

These in-depth simulations are necessitated by the long lead times and large amounts 
of money involved in constructing a supercollider detector. They have been made possible 
by the convergence of a reliable theoretical framework (the Standard Model) with rapid 
increases in available computing power. The simulations are an essential component of 
the design process, providing feedback on the design's strengths and weaknesses. In this 
way, the simulations have played and continue to play an important role in the design of 
detectors for the SSC and LHC. 

This article reviews the simulations carried out by the GEM (Gammas, Electrons, and 
Muons) Collaboration for its detector^^ at the SSC. Section 2 describes the tools available 
for simulation of interactions and of the behavior of the resulting particles in a detector. 
Section 3 briefly reviews the design proposed for the GEM detector at the time of Technical 
Review. It is essential to have a detailed engineering design to obtain realistic results, and 
the details of results presented here are, strictly speaking, valid principally for this design. 
However, any well-developed design is likely to have comparable strengths and problems, 
so the methods used for GEM and the results obtained should be more generally useful. 
Section 4 describes the tools developed for the simulation of the GEM detector. The prin- 
cipal development was a fast all-detector parameterization, called gemf ast, that allowed 
simulation of very large numbers of events in a manageable amount of time. This was sup- 
plemented by the most detailed possible simulations of the muon system and calorimetry in 
those cases that the fast simulation was insufficient to model important details and where 
the importance of the physics to GEM's design was believed to warrant the considerable 
extra effort. In section 5, the bulk of this article, we describe the results obtained for 
a variety of representative processes. These were chosen to impose stringent tests on all 
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areas of the GEM detector's design as well as for their intrinsic physics interest. Finally, 
Section 6 summarizes the results and conclusions. In publishing this document, it is our 
hope that others will see it as describing a useful model of detector simulation. We believe 
that that this model can be improved. We are also convinced of the validity of the simu- 
lation philosophy that culminated in gemf ast. This gave GEM the capacity to simulate 
rapidly a very wide variety of physics and background processes and to use the results of 
these simulations to provide input for crucial design improvements. Without gemf ast, we 
would have had to compromise either the breadth or the plausibility of our studies. 
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2. SIMULATION TOOLS IN HIGH ENERGY PHYSICS 

Since the Standard Model is a gauge theory, the couphng of any new physics to known 
particles is determined by the SU{3) (8) SU{2) U{1) representations of the new particles. 
Furthermore, since the SU (3) coupling becomes small at large momentum transfer and the 
SU{2) and U{1) couplings are intrinsically small, the cross sections for the production of 
new particles can be calculated in perturbation theory provided that the new physics does 
not involve nonperturbative dynamics. Even if the new dynamics is nonperturbative, its 
coupling to ordinary matter remains perturbative, and useful predictions can still be made, 
e.g., for strong WW scattering.^^ All the Standard Model background cross sections can 
likewise be calculated using perturbation theory. 

For testing a detector design, it is not sufficient to calculate cross sections; it is necessary 
to generate complete events; event generators are briefly reviewed in Section 2.1 below. It 
is also necessary to simulate the interactions of each of the particles from these events with 
the detector. One of the most generally useful programs for doing this is GEANT,^^ which 
contains tools for describing the geometry of the detector and for propagating particles 
through it taking into account all the possible interactions. CALOR89,^^ is less general 
but contains a better description of low-energy neutrons and so is better suited to precision 
calculations of the response of hadron calorimeters. These simulation tools are reviewed 
in Section 2.2 below. 

While GEANT and CALOR89 are very powerful, they are too slow to simulate the 
large number of events needed to determine the backgrounds for new physics signatures 
in hadronic collisions. Hence it is necessary to construct fast, parameterized simulations 
based on GEANT or CALOR89 results for single particles and on test beam data. A 
sophisticated parameterized simulation of the GEM detector, gemfast, is the basis of 
most of the physics results herein. It is described in Section 4 below. 

2.1. Review of Event Generation 

HERWIG,^^ ISAJET,^'^ and PYTHIA^^ are all programs widely used in particle physics 
to generate complete events for a variety of processes. While these programs differ in many 
important details, they share a common general approach, which is outlined here. 

The first step in producing a complete event is to generate an initial hard scattering 
process using its perturbative QCD or electroweak cross section. Because the QCD effective 
coupling q;s((5^) becomes small at large momentum transfer Q^, the cross section a for 
any hard process involving quarks and gluons can be calculated in perturbation theory. 
Then the pp cross section is given by the convolution of a with the parton distributions 

cr = y dxidx2 d- (s, Q^) S (^xiX2 - ^ fi {xi, Q^) fj {x2, Q^) ■ (2.1) 

Beyond the leading order in as, & in general contains infrared and coUinear singulari- 
ties. For inclusive cross sections, the infrared singularities cancel between real and virtual 



7 



graphs, and the coUinear singularities can be absorbed consistently into the parton distri- 
butions fi{XjQ'^)j producing their dependence. To lowest order, however, a is free of 
any singularities, so Eq. (2.1) can be used to generate an exclusive process using standard 
Monte Carlo techniques. That is, the program selects the four-vectors for the initial hard 
scattering. 

While Eq. (2.1) is a good approximation for sufficiently inclusive cross sections, it is a 
poor approximation for the structure of events. The hard scattering process represented by 
a involves external quarks and gluons at least in the initial state and generally in the final 
state as well. These are massless and carry color charge, so they necessarily radiate gluons. 
Since the coupling in QCD is dimensionless, this radiation occurs on all scales from 
that of the initial hard scattering down to the confinement scale, 0{1 GeV^). In particular, 
a massless quark or gluon can radiate a coUinear massless gluon while conserving energy 
and momentum, leading to a so-called mass or coUinear singularity. Because of these 
singularities, the higher order corrections are of order as{Q^)logQ^, and they cannot be 
neglected. 

Fortunately, the coUinear singularities take a simple form. Each additional radiation 
in which an external parton i produces j and k introduces into the cross section a factor 

:^P^^Jk{z) (2.2) 

where z = pj/pi is the momentum fraction carried by j and Pi^jk is one of the well-known 
Altarelli-Parisi functions. The dominant contribution comes from p'jj.<^Pi<^ Q^, so 
each radiation is approximately independent. Because of this and because the cross sections 
and not the amplitudes factorize, the quantum mechanical radiation in the coUinear limit 
can be treated as a classical Markov branching process. The radiation of soft gluons also 
factorizes, just as it does for QED. This makes it possible to write a Monte Carlo algorithm, 
known as the branching approximation, that correctly generates all the as{Q'^) logQ^ 
corrections to all orders in QCD perturbation theory. Each branching leads to partons with 
smaller pf. Since the Monte Carlo algorithm uses exact, non-coUinear kinematics, it also 
generates a fair approximation to non-coUinear radiation, i.e. to multijet final states. 

For sufficiently small pf perturbation theory becomes invalid, and one must resort to 
a model to convert partons into hadrons. Both ISAJET and PYTHIA use fragmentation 
models based on low-energy data for this. HERWIG carries perturbation theory down to 
pf = 0{1 GeV^) and uses mainly two-body phase space for the final hadronization. It is 
also necessary to model the hadronization of the spectator partons, which carry most of 
the total energy; there is no reliable theory for this. Fortunately, while the hadronization 
model is important for the details, perturbation theory controls the main features. 

An alternative to the branching approximation is to use explicit higher-order matrix 
elements. This approach produces a better approximation to cross sections for events 
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with several distinct jets, but it is difficult to combine with parton showers to produce 
realistic events without double counting. VECBOS^^ contains the exact tree- level cross 
sections for W + njets with n < 4, and has played an important role in searches for the 
top quark. PAPAGENO"^^ contains a wide range of processes and has been used in some 
of the analyses presented here. 

2.2. Review of GEANT and CALOR89 

GEANT is a system of detector-description and simulation tools written largely at 
GERN. The first version, written in 1974, was a bare framework which emphasized track- 
ing of a few particles per event through a relatively simple detector. GEANT has been 
continuously developed since then to handle much more complex detector systems, many 
hundreds of particles per interaction, and greater detail in the physics simulated. Be- 
cause of its flexibility, good support on many hardware platforms, and years of vcriflcation 
against data, GEANT has become a de facto standard within high energy physics and is 
used very extensively. 

GEANT performs the following functions: 

• Define particle and material parameters. GEANT has internal databases of 
particle masses and lifetimes, properties of common materials, and interaction cross 
sections. The user is able to define additional materials, and create mixtures or com- 
pounds. 

• Geometry definition. The user must create a geometrical model of the detector. 
GEANT will then handle the tracking of particles through this model and their inter- 
actions within it. GEANT 3.15 provides a palette of 14 standard shapes for the user 
which can be sized, positioned and rotated as desired. The detector model is a hier- 
archical tree in which each volume is contained within a 'mother', and each 'mother' 
contains one or more 'daughters.' 

• Compute physics processes. GEANT accepts lists of particles and vertices pro- 
duced by event generators. It contains routines to handle the decays of unstable 
particles, and the interactions of particles with the detector. The physics package 
computes the probability of occurrence of a process by sampling the total cross sec- 
tion and generates the final state after interaction. It also handles continuous or 
quasi-continuous processes such as energy loss. The processes currently implemented 
for electrons and photons are: multiple Coulomb scattering, ionization and delta ray 
production, bremsstrahlung, positron annihilation, e+e^ pair production, Compton 
scattering, photoelectric effect, photofission of heavy elements, and Rayleigh scatter- 
ing. For muons, decays in flight, multiple scattering, ionization and delta ray produc- 
tion, bremsstrahlung, direct pair production and nuclear interactions are implemented. 
For hadrons, decays in flight, multiple scattering, ionization and delta ray produc- 
tion, and hadronic interactions are included. The hadronic interactions are handled 
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by parameterizations of experimental data incorporated in the packages GHEISHA, 
HADRIN/NUCRIN33 and FLUKA. ^4 

• Track particles through the detector. Particles from the primary interaction, and 
those generated within the detector, are tracked through the detector model. Tracking 
proceeds in steps; during each step, the particle may decay or interact; otherwise, it 
loses energy through ionization and continues to the next step. Particles whose kinetic 
energy falls below a cutoff value are stopped and their remaining energy dumped. 

• Simulation of detector response. The detector response is modeled in two stages. 
Firstly the particles are tracked through the sensitive detector elements; then, after the 
whole event has been tracked, digitization routines (provided by the user) are called 
to model the response of the detector to the particles. This two-step process enables 
nonlinear effects (such as detectors able to respond only to the passage of one particle) 
to be properly modeled. 

• Drawing. GEANT includes a drawing package to produce graphical output showing 
the detector model and the trajectories of particles through it. Tools are provided for 
the user to add markers for the locations of energy loss, hit detector elements, etc. The 
drawing package can also produce graphical representations of the hierarchical tree of 
volumes or of the internal data structures of the program. 

Although GEANT is indispensable for general detector simulation, we have also made 
use of the CALOR89^^ package for calorimeter simulations. We have used CALOR89 both 
directly, and through an interface with GEANT. CALOR89 provides a more detailed low- 
energy treatment of hadronic processes, especially neutron transport, which is of interest 
in the prediction of hadronic calorimeter resolutions. It uses HETC with FLUKA34 for 
general hadronic transport, MORSE"^*^ for neutrons with energy below 20 MeV, and EGS^^ 
for electromagnetic particles. CALOR89 contains its own geometry, tracking and physics 
packages similar in concept to those of GEANT. 
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3. THE GEM DETECTOR 

The GEM collaboration was formed in 1991 to develop a major detector for the SSC. 
The physics objectives of GEM reflected the motivation for the SSC itself: to study high 
transverse momentum physics — exemplified by the search for Higgs bosons — and to 
search for new physics beyond the Standard Model. The GEM detector design was pre- 
sented in a Technical Design Report (TDR),^^ and the simulations described in this paper 
are based on that design (See Figure 1). The detector was designed as a general purpose 
device emphasizing clean identification and high resolution measurement of gammas, elec- 
trons, and muons, hence the name, GEM. High priority also was given in the design to 
accurate measurement of jets and missing transverse energy. The architecture and the de- 
tector technologies chosen were intended to provide good performance even at the highest 
luminosities that were contemplated for the SSC. The design of the GEM detector was 
based on the following principles: 

• Very precise electromagnetic calorimetry without a magnet coil in front of it. Together 
with the redundant pointing capabilities of the EM calorimeter and central tracker, 
this provides the best measurements of photon and electron energies and allows the 
reconstruction of the masses of narrow states with very high resolution. 

• A precise muon spectrometer in a large superconducting solenoidal magnet, allowing 
measurement of the momenta of high energy muons over a large solid angle with a 
minimum of multiple scattering. The muon system is situated in a quiet environment, 
shielded by the thick calorimeter. Muon measurement does not rely on central tracking, 
thereby providing a robust stand-alone system for operation at the highest luminosities. 

• Hermetic hadronic calorimetry for the measurement of jets and the reconstruction of 
missing energy. 

• Central tracking in a magnetic field with sufficiently low occupancy to operate reliably 
at the highest luminosities anticipated at the SSC (10"^^ cm"^ s~^). The central tracker 
radius is 0.9 m, allowing for a compact calorimeter and a large muon tracking volume. 
It was assumed in simulations that the inner silicon tracking layers were removed at 
£ = lO^^cm-^s-^ 

The performance specifications of the detector are given in Table 1. The reach of the 
GEM design for physics signals that were of of major interest to the SSC is summarized 
in Figure 2. 

3.1. Detector Design 

The detailed optimization of the design was determined by the physics requirements, 
the practical ability to meet the necessary performance specifications, and cost/schedule 
constraints. Attention was also paid to detector integration issues such as radiation shield- 
ing and the interaction between the beamline and the inner radius of the detector. It was 
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Scanned detector figure is too large for electronic distribution. 



FIG. 1: Perspective view of the GEM detector. 

also intended that GEM should be complementary to the SDC detector, which used a large 
tracker; GEM emphasized calorimetric measurement of gammas and electrons, a precise 
and robust muon system, and unique capabilities at high luminosity. 

Detailed descriptions of the detector subsystems, including technical features, imple- 
mentation and integration issues, and studies of expected performance are given in the 
Technical Design Report. We describe briefiy some of the key features of each major 
subsystem below. 

3.1.1. Magnet 

A very large superconducting solenoid surrounded the detector elements, with a nom- 
inal field of 0.8 T, an inner diameter of 18 m, and a length of 31 m. The magnet provided 
a nearly uniform axial field in the region of the central tracker, allowing measurements of 
the momenta and signs of charged particles and helping with electron identification. In the 
volume of detector outside the calorimeters, the magnet provided a 0.8 T field for muon 
momentum measurements. The coil was in two halves, one on each side of the central de- 
tector support (which carried the calorimeter and central tracker). In the forward region, 
conical iron poles were employed to shape the field. These introduced a radial component 
to the forward field, improving the resolution of the muon system in the forward direction. 
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Table 1: Top-level specifications for the GEM detector. 



Magnet 




Central field 


0.8 T 


Inner diameter 


18 m 


Length 


31 m 


Muon system 




Coverage 


0.1 < \r]\ < 2.5 



Apt/pt at 
Apt/pt at 



0, pt = 
2.5, PT 



500 GeV 
= 500 GeV 



Charge separation {r] = 0) 



5% 
12% 

p < 6.5 TeV at 95% C.L. 



Electromagnetic calorimeter 
Coverage 
Energy resolution 
Position resolution 
Pointing resolution 



\v\ < 3 

6-8%/v^ © 0.4% 
4.4mm/v^ 

40-50 mrad/ \/E + 0.5 mrad 



Hadronic calorimeter 
Coverage 
Jet resolution 



|?7| < 5.5 
60%/ © 4% 



Tracker 
Coverage 

Charge separation at 95% C.L. (77 = 0) 

Momentum resolution 

at high momenta (measurement limited) 

at low momenta (multiple scattering limited) 



\ri\ < 2.5 
p < 600 GeV 

Ap/p'^ = 1.2 X 10^3 GeV~^ 
Ap/p = 3.5% 



3.1.2. Muon System 

Muons were identified by their penetration through the calorimeter. Muon momentum 
was measured using the sagitta determined by three superlayers of chambers between the 
calorimeter and the magnet; see Fig. 3. Cathode strip chambers (CSC's) were chosen, based 
on the need to obtain the required spatial resolution, determination of the z coordinate, 
triggering, and occupancy, and the desire to use a single technology in both the endcaps 
and barrel. The resolution in the sagitta measurement varies as BL^, where B is the 
magnetic field strength and L is the lever arm of the measurement. Through the use of 
a large lever arm and a moderate field, a precise momentum measurement was obtained: 
Apt/pt = 5% (12%) at t] = (2.5) for pT = 500 GeV. 

The momentum resolution was limited at high momenta by the spatial measurement 
errors (both inherent and due to misalignment), and at low momenta by the multiple 
scattering in the middle superlayer of chambers and energy-loss fluctuations in the cal- 
orimeter. To be able to sign-select muons at very high momentum (e.g., from Z' decay 
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FIG. 2: The physics reach of GEM. 

at the highest mass, ~ 8-10 TeV, that was accessible at the SSC) demanded single layer 
resolutions of 75 fim. and alignment between superlayers of 25 /xm. To avoid degrading the 
low-momentum muons which contribute to the signal for ZZ* the middle 

muon chamber superlayer must be less than 10% of a radiation length in thickness. 

In order to be sufficiently robust to operate at the highest luminosities {C = 
10"^^ cm^^ s^^) attainable at the SSC, the muon system was shielded by a thick (> llA at 
r] = 0, increasing in the forward direction), nearly hermetic, calorimeter. The calorimeter 
thickness was chosen to keep the rate from punchthrough hadrons significantly below that 
from in-flight decay muons. In addition, a geometry for the forward calorimeter was chosen 
that minimizes the background. 

A feature of the muon system was the use of a 0.2 m open space between the calo- 
rimeter and the first muon superlayer. This clear space allowed charged particles from 
electromagnetic showers initiated by high-momentum muons to be bent out of the muon 
path, leading to increased reconstruction efficiency. 

3.1.3. Calorimeter 

The general layout of the GEM calorimeter is shown in Fig. 4. In order to identify and 
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FIG. 3: Quarter view of the GEM detector showing the muon system, including 
shielding. 

measure electrons and photons, a precision noble-liquid electromagnetic calorimeter was 
used, employing krypton in the barrel and argon in the cndcap. This technology was chosen 
because of its ability to achieve the required resolution, longitudinal segmentation and 
pointing ability, its intrinsic radiation resistance, its ease of calibration, and the extensive 
experience that has been acquired with large liquid-argon systems. An accordion geometry 
was used to provide good hermeticity and allow for faster readout than parallel-plate 
calorimetry because of lower inductance and capacitance. 

Behind the EM calorimeter was a hermetic hadronic calorimeter used both to measure 
jet energies and, in conjunction with the forward calorimeter, to measure IpT- The barrel 
hadron calorimeter was a hybrid system that performed the energy measurement primarily 
in a noble-liquid calorimeter (in the first ~ 6A), followed by a relatively inexpensive cop- 
per/scintillator calorimeter which provided energy measurement for late-developing show- 
ers and the necessary shielding for the muon system. In the end cap an all-liquid design 
was used. 

In the forward region, additional calorimeters measured high-momentum particles near 
the beam pipe. Together with the barrel and endcap calorimeters, these forward calorime- 
ters determined the sensitivity to ^t- In order to detect new particles such as massive 
gluinos and squarks, the forward calorimeter must cover the region |?7| ^ 5.5, be sufficiently 
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FIG. 4: Quarter view of the GEM detector showing the calorimeters. 

dense to contain hadronic showers, be sufiiciently fast to cope with the high particle flux 
in this region, and be radiation hard. The design had a first section consisting of a liquid- 
argon tube calorimeter, followed by a liquid-scintillator-capillary and tungsten calorimeter. 
The calorimeter was optimized to obtain good spatial information in the first section and 
adequate shower containment in the second section. 

3.1.4. Central Tracker 

The central tracker was 1.8 m in diameter by 3.5 m long, surrounding the interaction 
point. The tracker size was determined by the needs to obtain sufficient tt'^ rejection by 
shower shape analysis, to minimize the calorimeter cost, to optimize the tracker resolution, 
and to preserve sign-selection ability to high momenta. The layout of the tracker system 
is shown in Fig. 5. 

The tracker supported the goal of measuring gammas, electrons, and muons at high 
pr, at high luminosity, up to £ = 10^^cm~^s~^. The tracker provided good separation of 
gammas and electrons by finding a charged track, and measured the electron sign up to 
600 GeV. The former is essential to the search for — > 77 and to background rejection 
in Z' — > e+e~; the latter, for the gluino search using the signature of same-sign leptons. 
Another important role for the tracker was to measure the position of the primary vertex. 
Physics involving b, t, and r decays requires full pattern recognition capability, including 
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FIG. 5: The GEM central tracker. 

secondary vertex finding and tracking at low momenta. 

The design used two technologies. For the inner section of the tracker, silicon mi- 
crostrips were used, giving the required fine segmentation and radiation resistance. For 
the outer section, interpolating pad chambers (IPC's) were chosen for their low occupancy, 
their ability to provide near three-dimensional space points, their high-luminosity capabil- 
ity, and their demonstrated operational resolution of 50 /im. 

3.1.5. Electronics/Data Acquisition 

Triggering and data acquisition in GEM followed a three-level strategy to provide a 
system without deadtime that provided as much information as possible at each trigger 
level. It was designed for luminosities up to 10^^cm~^s~^, with provision for operation 
at higher luminosities with modest upgrades. The architecture consisted of a synchronous 
and pipelined Level 1, an asynchronous Level 2 (possibly with special purpose hardware), 
and a Level 3 processor ranch. Full granularity data was available at Levels 2 and 3. Level 

1 was designed to handle up to 60 MHz input rate, with an output rate of 10 kHz. Level 

2 was designed to handle an average input rate up to 100 kHz, with an output rate of 300 
Hz. Finally, Level 3 accepted 3 kHz, with an output rate of 100 Hz. 
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4. MODELING THE GEM DETECTOR 

In principle, the full detector simulation programs GEANT and CALOR89 could be 
used to compute GEM's response to the signal and background for any process. However, 
these programs are very slow for complex events at high energy. It is impractical to use 
them to simulate the more than 10^ events often needed to determine a rare signal's back- 
ground arising from a combination of relatively likely processes. Consequently, two kinds 
of simulation were done for GEM. Detailed simulations, based on GEANT, of each of the 
individual detector systems have been performed for single particles or for limited numbers 
of complete events. The results of these detailed studies have been parameterized and 
incorporated in gemf ast, the fast simulation program for GEM that was used for determin- 
ing the performance of the detector for physics processes. For example, parameterizations 
were made of energy and momentum resolutions and of muon reconstruction efficiency in 
the presence of other particles in the event of interest. Most of the results presented in this 
article are based on gemf ast. Where necessary, hybrid simulations of gemf ast and full 
GEANT have been used. For example, in the study of H'^ — > 77, detailed electromagnetic 
shower shape studies for real photons and and jets faking photons were carried out with 
GEANT. Apart from examples such as this, gemf ast describes the performance of GEM 
quite accurately. 

The key to a fast detector simulation is to use a very simple geometry and to param- 
eterize the response of each detector component in a simple way. The geometry used in 
gemf ast is a set of concentric cylinders, one each for the central tracker (CT), electromag- 
netic calorimeter (EC), hadronic calorimeter (HC), scintillator calorimeter (SC), forward 
calorimeter (FC), and muon system (MU). The geometry is shown in Fig. 6. The radia- 
tion and absorption lengths in the cylinders representing the calorimeters are varied with 
r] so as to match the true design. Particles are tracked through each successive volume on 
straight lines for neutral particles or on helixes in a uniform magnetic field, B = 0.8 T, 
for charged ones. Since the MU resolution is parameterized, particles need not be tracked 
through the nonuniform field in the muon system. 

Once a particle enters a given detector system, its energy resolution, angular resolu- 
tion, and detection efficiency are calculated using parameterizations of the single particle 
response. This simple, single-particle approach is not adequate for the central tracker 
reconstruction efficiency, which is sensitive to the presence of other tracks in the same 
event and to pileup, and which therefore has been investigated separately. Development of 
electromagnetic and hadronic showers in transverse and longitudinal directions, including 
fiuctuations, is modeled and the individual tower energies are calculated. Unstable parti- 
cles are allowed to decay anywhere in the detector using code adapted from the GEANT 
package. The Level 1 trigger response is also simulated. 

4.1. Event Generation 

The first step in simulating a process is to generate the events of interest. This is 
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FIG. 6. Geometry used in gemf ast. The detector is approximated by a set of concentric 
cylinders. Particles are tracked through the central tracker and calorimeters assuming a 
uniform magnetic field. 

handled by gemgen,^^ which incorporates ISAJET,^'^ PYTHIA,^^ and a single particle 
gun. The single particle gun generates a single particle at a given pT, r], and or a single 
quark or gluon fragmented with PYTHIA. 

The event generators are interfaced to the detector simulation in a way that allows 
adding a signal event and a Poisson-distributed number of minimum bias events in the 
same bunch crossing. The same or different generators can be used for each sample. The 
vertex position of each event is generated according to the expected width, cr^ = 5 cm. The 
same approach could be used to describe pileup from out-of-time bunches. However, this 
is impractical because of the computing time required to generate minimum bias events 
over the tails on the sensitive time of the detector. Instead, out-of-time pileup events have 
been generated separately, the detector response to them has been parameterized, and the 
resulting noise or inefficiency taken into account for each of the detector systems. 

4.2. Central Tracker 

The single particle momentum and vertex position resolutions of the central tracker 
were calculated using a full GEANT simulation including the effects of the magnetic field, 
detector geometry, chamber positions and resolutions, material distribution including sup- 
ports and cables, geometrical acceptance, detector efficiency, and distribution of the inter- 
action vertex. The effect of out-of-time pileup was included as an additional inefficiency 
due to the detector deadtime. Both the silicon strips and the IPCs are used at a luminosity 
of 10'^^ cm~^s~^; the silicon detector was assumed to be removed for 10^^ cm~^s~^. Figures 
7 and 8 show the resulting parameterization of the mean muon or pion resolution vs. pT 
and r] used in gemf ast for 10^^ cm^^s"-*^ and 10^^cm~^s~^, respectively, assuming a vertex 
constraint. The impact parameter resolution for tracks is also parameterized, allowing 
modeling of vertex reconstruction on an event-by-event basis. 

The reconstruction efficiency for isolated tracks exceeds 97% over the whole range of 
pseudorapidity \r]\ < 2.5. This is based on requiring at least ten good hits on a high-py 
track. The number of hits is calculated in gemf ast using the positions of the chambers and 
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FIG. 7. Parameterization of vertex-constrained central tracker resolution in gemf ast for 
muons or charged hadrons vs. pT and r] ior C = 10^^ cm~^s~^. Both the silicon tracker 
and the interpolating pad chambers are used. 
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FIG. 8. Parameterization of vertex-constrained central tracker resolution in gemf ast for 
muons or charged hadrons vs. pT and rj for C = 10^^cm~^s~^. The silicon tracker is 
assumed to be removed. 

the actual origin of the track. It is difficult to parameterize the reconstruction efficiency 
for non-isolated tracks, so the simulations described here use only tracks which have pt > 
1 GeV and which are isolated at the generator level in a Arj x Acf) region corresponding to 
three pads in the IPCs. 

The momentum resolution for electrons has been treated separately, taking into ac- 
count the emission of bremsstrahlung photons caused by the material in the tracker. (The 



20 



30 

1 20 

cn 

c 

> 
LU 

10 



-5 5 -1.5 1.5 

d(1/Pmeasured " ^^P) >< P d(1/P,,,3,,,, - 1/P) X P 

TIP-03946 

FIG. 9. Calculated momentum resolution function for electrons in the central tracker at 
Pt = 525 GeV for < < 0.2 (left plot), and pt = lOGeV for 2.2 < r] < 2.4 (right plot). 
The histogram shows the GEANT simulation, and the smooth curve is the parameteriza- 
tion used in gemf ast. 

photons are assumed to be emitted nearly parallel to the electron and to hit the same 
calorimeter cell, so that the calorimeter energy resolution is not degraded.) Figure 9 shows 
the electron momentum resolutions at a particular value of pt and i] from the full GEANT 
simulation and from the corresponding parameterization used in gemf ast. The parameter- 
ization fits the GEANT data well, including the bremsstrahlung tails. These are important 
for electron sign determination, which can be done up to pT ~ 600 GeV. 

Material in the central tracker is parameterized as a function of i] and (p after each 
layer of the silicon or IPCs, and this is used to convert photons at appropriate space points. 
Secondary e'^e~ pairs are generated using code adapted from GEANT. 




4.3. Calorimeter 

The parameterized response of the central calorimeter, which covers \r]\ < 3, includes 
shower profiles and energy resolutions. Longitudinal and transverse electromagnetic and 
hadronic shower profiles are generated using GFLASH 1.3, which was originally devel- 
oped to describe the HI liquid argon calorimeter.^^ It has been modified to work with 
the gemf ast geometry. GFLASH incorporates correlated fluctuations of shower profile 
parameters, hadronic shower fluctuations into early tt'^'s, transverse proflle variations with 
depth, and shower development along the true direction of incidence as determined by 
tracking through the central tracker region. It gives a good description of the shapes 
of both electromagnetic and hadronic showers in uniform regions of the calorimeter, as 
illustrated in Figure 10. GFLASH is used to distribute shower energies among calorim- 
eter towers. The electromagnetic, liquid-hadronic, and scintillator tail-catcher towers are 
modeled separately, but further longitudinal segmentation is neglected for faster execution. 
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FIG. 10. Ratio of electromagnetic to total energy for 100 GeV pions. The histogram shows 
the result of a full GEANT simulation, and the points show the GFLASH-based parame- 
terization used in gemf ast. 

The segmentation of the calorimeter is varied realistically with rj. The simulated 
endcap segmentation is in rj-cf) rather than in x-y, but this should have no effect on physics 
performance. The EC and HC segmentations are Ar] x A(j) = 0.026 x 0.026 and 0.08 x 0.08 
in the middle of the barrel. Since the segmentation is approximately constant in units of 
radiation or absorption lengths, it is about a factor of six coarser at the small-angle edge 
of the endcap than in the middle of the barrel. 

The energy resolution of the electromagnetic calorimeter has been calculated using a 
GEANT simulation including an extremely detailed geometry and very low cutoffs. A sim- 
ilar simulation gave good agreement with test beam data for the non-projective accordion. 
The resolution has been parameterized in the form 

AE a{r]) 



(4.1) 



E y/E 

with the parameters tuned to give the correct resolution as a function of E and ry for the 
5x5 sum of towers used to obtain optimal resolution for isolated electrons and photons. 
Typically, the stochastic term (a) is about 6% in the barrel and about 8% in the endcap, 
and the constant 6-term is about 0.4%. Pileup and noise are not included in the energy 
resolution but are added separately. 

The pointing resolution of the electromagnetic calorimeter for photons was calculated 
using the same detailed GEANT simulation. It is parameterized as 



A0 



sfE 



(4.2) 



for low to moderate energies, where is about 40 mrad in the barrel and 50 mrad in the 
endcap and 651(77) is about 0.5 mrad. At the highest energies, for which the shower is no 
longer fully contained in the EM calorimeter, 2 mrad. 
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FIG. 11. The histogram shows the distribution of pileup, including out-of-time events, in 
gemf ast for a 0.16 x 0.16 electromagnetic trigger tower. 

Thermal noise is added to each tower depending on rj. In-time pileup events are added 
explicitly. To determine the effects of out-of-time pileup, a Poisson-distributed number 
of minimum bias events is generated for each of the 50 preceding and 20 following bunch 
crossings and simulated with gemf ast. The calorimeter response for each bunch crossing 
is weighted with a response function that takes into account the intrinsic calorimeter 
response and the shaping circuit in the electronics, using shaping times of 40 ns for the 
electromagnetic calorimeter and 100 ns for the hadronic calorimeter. The sum for each cell 
is calculated to provide a snapshot of the response to pileup. This ensures that all the 
longitudinal and transverse correlations among cells, caused either by individual shower 
shapes or by jets, are preserved. One of these snapshots is then superimposed on the 
response from the signal and in-time pileup events. As can be seen from Fig. 11, this 
approach gives larger fluctuations and a smaller half-width than the equivalent Gaussian 
noise. For the analyses described here a sample of 10^ such snapshots has typically been 
used. 

The GEM calorimeter design has an intrinsic e/h ^ 1, so reconstructing a jet energy 
as a simple sum of the observed energies in the various layers would give a large constant 
term in the resolution. This effect can be reduced by using an iterative weighting proce- 
dure. ^^''^'^ To save execution time this weighting is not implemented in gemf ast. Instead, 
the sampling and constant terms in the single hadron resolution have been tuned so that 
the jet energy resolution from gemf ast matches that from the detailed GEANT simulations 
including the weighting. The resulting resolution for isolated jets is 



AE 



E 



0.6 



jet 



®0.04. 



(4.3) 



For many cases the effects of the clustering algorithm used to define jets are comparable 
to those of energy resolution. 
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FIG. 12. Calorimetric Et resolutions for hadrons vs. rj ioi E = 10 GeV (solid line), 50 GeV 
(dashed line), 200 GeV (dotted line), and 500 GeV (dash-dotted Une) used in gemfast for 
the calculation of ^t- Most of the fluctuations are caused by limited statistics in the 
GEANT simulation. 

The forward calorimeter covers 3 < |?7| < 5.5, with full measurement capability to 
~ 5.0. It has been used in the physics studies described here only to determine the 
missing energy, ^j-. For this, detailed simulation of the response of individual cells is not 
needed. Rather, the energy and direction of each particle hitting the forward calorimeter 
is smeared according to a parameterization derived from a GEANT simulation. The sim- 
ulation includes all the effects of dead material and shower spreading across calorimeter 
boundaries. The resulting Et resolution as parameterized in gemfast is shown in Fig. 
12. The statistics in this simulation were not sufficient to study potential nongaussian 
tails. These are modeled in gemfast by adding a second Gaussian with a small amplitude 
and a larger width. Test beam data for single hadrons from 50 to 100 GeV in the D0 
liquid argon calorimeter show a tail composed of roughly 1% of the events with a standard 
deviation two to three times larger than the Gaussian calorimeter resolution. This tail 
is slightly larger than that seen for 1 TeV jets in GEANT studies of jet resolution using 
energy-dependent weighting. A similar tail has been assumed for GEM. 

4.4. Muon System 

The muon momentum resolution has been calculated by a full GEANT simulation 
including a detailed model of the detector and its support structures. The calculation 
includes chamber resolutions and alignment errors, the calculated shape of the magnetic 
field, the number of CSC planes in the measurement, and multiple scattering from the 
chambers and their supports. It is described at length elsewhere. The resulting param- 
eterization of the resolution is shown in Fig. 13. The jumps in the curves come from the 
transitions between various sets of chambers in the barrel and endcap. The flattening in 
the resolution at 77 ~ 2.5 is the effect of the forward field shapers. 
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FIG. 13. Parameterization of muon resolution vs. r] for pT = 100 GeV (squares), 500 GeV 
(circles), and 1 TeV (triangles) in gemfast. 



1.25 




FIG. 14. Parameterization of muon acceptance vs. r/ in gemfast. Because of binning of the 
GEANT data, the gaps in the barrel are wider than in reality but have nonzero acceptance. 
The average acceptance is correct. 

The geometrical acceptance for muons has been calculated similarly, requiring that the 
muon pass through at least three chambers in each superlayer. (This is overly conserva- 
tive,^^ since fairly good muon resolution can be obtained from two superlayers at low pj^ 
or from two superlayers plus the vertex at high pT-) The three-super layer acceptance is 
shown in Fig. 14 and is essentially independent of muon energy. The coverage of the region 
\ri\ < 2.46 is 83%. The losses due to the spoiling of hits by 6 rays and other electromagnetic 
interactions have also been simulated by GEANT, parameterized as a function of muon 
energy, and included in gemfast. 

The muon energy loss in the calorimeter has been calculated, including all the processes 
modeled in GEANT, and has been parameterized with a simple analytic function. The full 
calculation and the parameterization are shown in Fig. 15. The lost energy is added to the 
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FIG. 15. Muon energy loss in calorimeter for = 500 GeV (solid curve) and gemfast 
parameterization (dots). Note that the tail at large energy loss is fit well by the parame- 
terization. 

appropriate cells of the calorimeter so that energy losses large compared to the noise can be 
reconstructed in gemfast. The best resolution is obtained by using the measured energy if 
it is greater than a factor k ~ 1.5 times the most probable value, and the truncated mean 
value otherwise. 



4.5. Trigger Simulation 

The Level 1 trigger efficiency is simulated in gemfast and included in the analyses 
described below. The Level 2 and Level 3 triggers are assumed to be included, to a first 
approximation, in the cuts made in the analysis. 

An initial set of trigger thresholds is shown in Table 2. The Ex triggers include an 
isolation cut on the surrounding towers and on the hadronic energy behind the cluster. The 
Jx triggers sum the electromagnetic and hadronic sections in towers reasonably matched to 
the size of a QCD jet. The Mx thresholds are nominal values at which the acceptance is 84% 
and include the coarse resolution of the trigger roads. The missing energy simulation uses 
the full sum of the calorimeter. While the highest thresholds in Table 2 have acceptable 
rates by themselves, the lower thresholds can only be used either with prescaling or in 
combinations. 

Twelve combinations of the primitives listed in Table 2 appear to be sufficient to select 
all the physics processes so far considered for GEM; they give a total Level 1 trigger 
rate less than the design goal of 100 kHz. In particular they provide triggers with good 
efficiency for low-mass Higgs bosons and the other (relatively) low-mass processes listed 
in Table 3. The most difficult trigger is that for inclusive — > 77 at the low end of the 
mass range. Some low-pT processes such as jets have such large cross sections that they 
must be prescaled. No attempt has been made to implement an efficient trigger on the 
inclusive-6 cross section to study 6-physics, although this might be of interest during the 
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Table 2. Primitives used in the GEM trigger. The highest threshold can be used alone; 
the rest can be used in conjunction with other triggers. 



Name 


Definition 


Thresholds 


Ex 


EM cluster with Et > x GeV 


Eg, EiQ, £50, Ego 




in 0.16 X 0.16, isolated in both 






EM and HAD calorimeters. 




Mx 


Muon with pT > x GeV. 


Mio, M20, M30, M40 


Jx 


Jet cluster with Et > x GeV 


J16, J50, Jso, J200 




in 0.48 X 0.48. 




$x 


Missing energy > x GeV 


-^505 -^100 




summing whole calorimeter. 





Table 3. Trigger efficiencies for a variety of processes using the triggers shown in Table 2 
with standard logical notation. The triggers listed are the principal ones among those 
studied, but the efficiency is for the sum of all combinations. The efficiency for H — > 
£~^£~£~^£~ at 140 GeV could be increased about 5% by adding an Eiq A Mio trigger. The 
last four processes are useful mainly for calibration. 



Process Mass Trigger Efficiency 

80 GeV 2EiQy E^q 78.7% 

iy^77 140 GeV 2EiqV E50 94.7% 

ttH^^^£X 80 GeV 2Ei6 94.4% 

£+£-£+£- 140 GeV 2Mio V 2£;i6 81.8% 

£+£-£+£- 400 GeV 2Mio V 2£;i6 99.8% 

H^£+£-jj 800 GeV 2MiQy 2Eiq 99.9% 

tt^£ubX 140 GeV E50 V M30 75.3% 

gg 500 GeV 3 Jso A ^100 99.9% 

W^eiy — EiqA$^q 15.8% 

W ^ ixv — Mio A -^50 48.7% 

Z-^ee — 2Ei6 80.3% 

Z fifi — 2Mio 86.9% 



low-luminosity initial period. Level 1 triggers on the higher mass processes which have 
been considered are generally easier than those listed in Table 3. 
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5. STUDIES OF PHYSICS PROCESSES 



5.1. Selection of Processes 

Wc have studied the performance of the GEM detector for a variety of physics processes 
chosen both because they probe the capabihties of important components of the detector 
and because of their intrinsic physics interest: 

• Section 5.2 presents an in-depth study of the search for the Higgs boson of the Stan- 
dard Model. The signals, backgrounds and discovery potential for M}j = 80 — SOOGeV 
are discussed. Depending on the Higgs mass, the modes studied are — > 77; 
H^{tt/W) £±77X; ^ ZZ^*) 4 charged leptons; ^ ZZ ^ e+e'Vu; 
and ^ ZZ ^ jetjet. 

• Flavor physics involving top-quarks is discussed in Section 5.3. We describe the mass 
measurement of a heavy top-quark in the standard decay mode t — > W'^b using two 
methods: t — > isolated plus non- isolated //~, and t — > 3 jets. We also discuss the 
discovery of a charged Higgs boson in the nonstandard decay mode t — > H^b, followed 

by H+ T+Ur. 

• Jet physics is studied in Section 5.4. We discuss the determination of the jet energy 
scale, using as a physics context the search for quark substructure in high-E'j- jets. 
Other jet studies are discussed in Sections 5.2 Z^Z^ £+^-jetjet) and 5.3 
{t^W+b^ 3 jets). 

• Supersymmetry is discussed in Section 5.5 as an example of missing transverse energy 
signatures for new physics. The distribution is calculated for GEM, including the 
effects of transition regions and dead material. The signature is studied for a range 
of gluino and squark masses. In addition, the like-sign dilepton signature for gluino 
production is investigated. 

• Section 5.6 is devoted to studies of high-mass-scale physics at ultrahigh luminosity (for 
the SSC), C ~ 10'^'* cni^^ s~^. These physics studies include precision investigations 
of the properties of a very massive Z' boson in its e+e" and n^/j,^ decay channels, 
and the character of quark-lepton substructure contact interactions via the process 

qq fi+fi-. 

For all these processes, the performance of the GEM design has been determined 
realistically. This is an important step in optimizing the design. Results are given below 
for a variety of integrated luminosities: 10 fb~^ and 30 fb~^, corresponding to one and three 
years of operation at the design luminosity of 10^^cm~^s~^; and 100 fb~^, corresponding 
to one year at a peak luminosity of 10^^cm~^s~^. 

As noted above, the upper limit on the top-quark mass was 95 GeV at the time most 
of the work described here was done. Except in Section 5.3, a top mass of 140 GeV was 
used in the simulations. If a top mass of 175 GeV were used, this would sometimes increase 
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signals (as in gluon fusion of a standard Higgs boson) and/or backgrounds and sometimes 
decrease them. Overall, we do not think that a significant quantitative change would result 
from using the higher mass. 

5.2. Standard Model Higgs Physics 

Here we describe searches for the minimal Standard Model Higgs boson with the GEM 
detector in the mass range between 80 and 800 GeV. This covers the interval between 
lower limits from LEP (which have excluded the mass range below ~ 60 GeV,^^ and should 
reach 80-90 GeV at LEP 200;^°and the triviality bounds discussed above.^^ The details of 
these Higgs boson search studies are presented in Ref. 51. Higgs bosons in extensions of 
the standard model such as the minimal supersymmetric model^^ generally have similar 
signatures, albeit with different cross sections. 

The dominant decay modes of the standard model Higgs boson are VF+VF" and ZZ 
for Mn > 2Mw and heavy fermion pairs for Mh < 2Mw- The latter, and the dominant 
decays of W^s and Z's into jets, all have large backgrounds, so it is necessary to rely on 
rare decays. The cleanest channel is — > ZZ/ZZ* — > {£ = e, which has 

four isolated high-pr leptons in the final state. This channel has an inadequate rate at the 
lowest and highest ends of the mass range. Hence, the search may be divided into three 
mass regions: 

• A Higgs boson of "intermediate" mass (80 GeV < Mh < 2Mz) will be searched for 
through its decays — > 77 and H ZZ* i^i~£'^£~. Both direct i!/" — > 77 and 
lepton-associated production {tt/W)H ^^£X will be used. 

• A heavy Higgs {2Mz < Mjj < 600 GeV) will be searched for through the channel 
H ^ ZZ ^ £+£-£+£-. 

• A very heavy Higgs {Mh 800 GeV) will also be searched for in £'^£~£'^£~ and in the 
channels H ^ ZZ ^ £+£-vv and H ^ ZZ ^ £^£~jj- 

The Monte Carlo event generator used in this study was PYTHIA 5.6 and JETSET 

7.3. ^^ The top quark mass was assumed to be 140 GeV. The study was carried out for the 
SSC design luminosity C = 10^^cm~2s~^ 

Table 4 summarizes the significance with which GEM could discover the Higgs through 
the channels listed above. This demonstrates the sensitivity of the GEM detector at greater 
than the 5a level for the whole range of Higgs masses considered. The signal can be seen 
with an integrated luminosity of 10fb~^, except for Mh < 100 GeV where 20 to 30fb~^ 
are needed. 

The most difficult part of the Higgs spectrum is the low end (between 80 and 100 GeV) 
and the region between the WW and ZZ thresholds, 170-180 GeV. GEM is sensitive to a 
Higgs with Mh below 100 GeV in both the two-photon and in the two-photon plus lepton 
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Table 4: Sensitivity of the GEM detector to standard model Higgs boson signals. 
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channels. This would allow the confirmation of a Higgs signal in two decay channels 
necessary for a credible discovery claim in this mass region. 

Over the entire Higgs mass range the dominant production mechanism at the SSC 
is gluon fusion, although vector-boson fusion is also important. The Higgs can also be 
produced associated with a vector boson or a top quark pair (ti) which provides an isolated 
lepton tag. 

In the standard model, a Higgs boson with mass Mh > 2Mz has a width which grows 
like M^, but it remains fairly narrow for Mh < 600 GeV. Figure 16 shows the lowest order 
Higgs production cross section multiplied by specific decay branching ratios for the five 
modes studied in this section. These cross sections were calculated using PYTHIA 5.6 with 
the EHLQ-1 parton distribution functions and a top quark mass of 140 GeV. Calculations 
of these cross sections using the more modern HMRS^^ ^nd CTEQ^^ distributions differ 
from EHLQ-1 by less than 20%.^^ Higher order QCD corrections are not included in 
Fig. 16. 

The physics signatures of Higgs decays all involve identification of isolated photons, 
electrons and muons. Since the Higgs cross sections are small and the potential back- 
grounds are large, this identification requires particular care. First, an isolation cut is 
applied, removing most of the jet background. Then, a detailed identification algorithm 
is used. The isolation cut for selecting electron and photon candidates was done using 
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Mh (GeV) 



FIG. 16: Lowest order Higgs cross sections into observable modes vs. mass, 
gemf ast and requiring 

^Et- E'^J^ < Ef^^"" + Ef^ , (5.1) 
R 

where the sum is over calorimeter towers in a cone of radius R. Generally, R = 0.45 was 
used for two-body final states and R — 0.30 — 0.35 for four bodies. The transverse energy 
Ej/' of the 7 or e candidate was found by summing the energy deposited in 5 x 5 cells 
in the electromagnetic calorimeter; i?™^^'^ is the mean transverse energy from pileup and 
noise; and E^^ is the isolation threshold imposed. 

The E^^ value was determined from the distribution of the thermal and pile-up noise 
as described in Section 4.3. A Gaussian fit to the total noise in an — cone of radius 
R — 0.45 gives a width of 3.4 GeV. This is reduced to 1.7 GeV if only those cells with 
\Et\ > 0.5 GeV are summed. Doing this shifts the E^^^'^ from 0.22 to 1.5 GeV, and reduces 
the signal loss from 13% to 7.6% for E^"^ = 4 GeV, the value used here. 

For photons, there must be no reconstructed charged track in the 5x5 cells. The 
central tracker can reject 95% of the electrons while keeping 96% of the photons. For 
electrons, there must be exactly one charged track in the central 3x3 cells. Also, the 
energy and shower position measured in the EM calorimeter must match the momentum 
and track position measured in the central tracker. The background to ilf — > 77 from 
misidentified Z — > e+e" is potentially large for Mjj Mz and is discussed below. 

A more complete understanding of the capability of GEM to identify isolated photon 
and electron candidates based on their shower shape requires a detailed GEANT simula- 
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tion. A sufficiently detailed geometry was available for only three 11 x 11 cell regions of 
the barrel. Therefore, the photons hitting the 5x5 cluster were rotated to the nearest 
of these regions, and the full simulation was run. Shower shape cuts (including lateral 
and longitudinal shower distributions), information from the narrow 9 strips of the first 
segment of the EM calorimeter, and energy in the hadronic calorimeter are combined in 
a likelihood function and used to select single electromagnetic showers. Typically this re- 
jects 75% of the remaining jet background while keeping 90% of the single photons. The 
probability of a QCD jet faking an isolated electron, 7?.(e/jet), was determined to be ap- 
proximately 10~^ by using detailed GEANT simulation. ^^'^^ For channels needing very 
good resolution, such as 77 and H ZZ* — > e^e~e^e^ , photons and electrons were 
required not to be in the region of degraded resolution between the barrel and the endcap, 
1.01 < |?7| < 1.16. The overall photon identification efficiency, including this geometrical 
loss and other cuts, is between 80% and 85%. The electron identification efficiency is 85% 
- 90%. For signals such as H ^ ZZ — > e+e^e+e^ for Mjj > 2Mz, one electron was 
allowed to hit this transition region. 

Muons were simulated using gemfast. They are relatively free of jet background and 
so require less stringent identification cuts. An isolation cut was made. 



where AE is the energy loss of a muon in the calorimeter, calculated as explained in 
Section 4.4. This eliminates muons from b and c decay, secondary decay muons in jets and 
punchthrough. The isolation cone radius was taken to be 0.35 for intermediate mass Higgs 
searches and 0.3 for heavier masses. The angle and momentum of the track measured in 
the muon system is then matched to that in the central tracker. The muon identification 
efficiency is around 80% for muons within \r]\ < 2.5, which includes factors of 85% from 
the geometrical acceptance and 95% from the muon track reconstruction and identification 
efficiency. 

Significance is used to indicate how well a Higgs signal can be identified in the presence 
of background. As in Gaussian statistics, a probability of 1 — 1.35 x 10~^ is expressed as 
a significance of 3(7, 1 — 2.85 x 10~^ is expressed as 5cr, and so on. In this section, a 5cr 
significance is generally regarded as the minimum for discovery. 

5.2.1. if ^ 77 Search for 80 GeV < Mr < 160 GeV 

Precision electromagnetic energy resolution is essential for the if — > 77 search because 
of the small production cross section (60 to 200 fb), the narrow decay width (3 to 100 MeV) 
of the Higgs boson between 80 and 160 GeV, and the large QCD 77 background. There 
is a potentially much larger background from QCD jets fragmenting into electromagnetic 
particles, which can be reduced using isolation cuts and a shower shape analysis made 
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possible by the fine segmentation of the calorimeter. In PYTHIA, jets can radiate prompt 
photons as an option. This is properly regarded as an approximation to part of the higher 
order QCD corrections to the 77 background. Since these corrections are explicitly known 
for both the Higgs signal and the qq 77 background, the prompt photon radiation has 
been turned off, and the PYTHIA signal and background cross sections have been scaled 
to the higher order results. 

Events were selected using the following cuts: 

1. \r]^\ < 2.5 and > 20GeV. 

2. Photon isolation cut with R = 0.45 and E^^ = 4GeV, and veto on 1.01 < |?7| < 1.16. 

3. Photon identification based on detailed shower shape. 

4. Electron rejection, as described below. 

5. I cos 0*1 < 0.7, where 9* is the polar angle of photon in the center of mass system of 
two photons. 

Cut 1 simply ensures that the photons are in the overall acceptance of the detector and are 
triggered on. Cuts 2 and 3 reduce the large potential backgrounds from misidentified QCD 
7-jet and jet-jet events well below the 77 continuum. Cut 4, which is only important for 
Mff f» Mz, removes the background from misidentified electrons. Finally, cut 5 reduces 
the real and fake 77 backgrounds. This leaves the iif — > 77 signal as a narrow bump in 
the Mjj distribution. 

The significance of a ilf — > 77 mass peak is directly related to the 77 mass resolution 
AM^^. This is given by 



where Ei and E2 are the energies of the two photons and 9 is the opening angle between 
them. It is clear from this equation that uncertainties in both energy and direction mea- 
surement would degrade the Higgs mass resolution. The energy resolution is determined by 
the EM calorimeter. The angular resolution is determined both by the position resolution 
in the EM calorimeter and by the precision of the vertex determination. GEM's calorime- 
ter position resolution. Ax — 4.4mm/\/E, has a negligible effect on the mass resolution. 
The vertex position for a single event is well determined by the central tracker. The only 
issue, therefore, is how well the correct vertex can be selected in the presence of an average 
of 1.6 additional minimum bias events at the standard SSC luminosity plus events from 
previous bunch crossings. 

Two approaches to determining the correct vertex are possible. The first is to use 
the difference in event topology resulting from the fact that Higgs production is a harder 
process than most of the minimum bias events and hence radiates more gluons. This leads 
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FIG. 17. I cos 6'* I distribution for an 80GeV Higgs, NLL direct photon background, 7-jet 
and jet-jet background. The cut used in this analysis is | cos^*| < 0.7. 

to a higher multiphcity and a higher average pt- Selecting the vertex with the highest 
pr-weighted charged multiplicity gives the correct Higgs vertex with 95% probability at 
the standard SSC luminosity. An independent approach is to use the pointing provided 
by the longitudinal segmentation of the calorimeter; using this pointing alone without any 
information from the central tracker degrades the mass resolution by only 20%. Using the 
pointing and then selecting the closest central tracker vertex gives a vertex within 5 mm 
87% of the time at standard luminosity. This degrades to about 65% at 3 x lO"^"^ cm~^s~^. 
A combination of these two methods improves the vertex-finding efficiency to 97%.^^ A 
vertex finding efficiency of 95% was assumed for this analysis. 

The overall M^^ resolution varies slowly from 0.66 GeV at M^^ = 80 GeV to 0.99 GeV 
at Myy = 160 GeV. Over the same range the Higgs natural width varies from 0.003 GeV 
to 0.097 GeV. 

The H — > 77 signal and the qq 77 and gg 77 backgrounds were generated 
with PYTHIA. The next-to-leading-log (NLL) corrections both to the signal^^ and to the 
qq — > 77 background^^ have been calculated. They enhance the signal cross section by 
a factor K Ki 1.5 and the background by a somewhat larger factor. The PYTHIA cross 
sections have been rescaled to these NLL results. 

After cuts 1 to 4, the cut on cos^* provides an additional rejection of both direct 
photon and jet backgrounds and thus improves the significance. Figure 17 shows the 
cos 9* distribution for an 80 GeV Higgs, NLL direct photon background and the 7-jet and 
jet-jet backgrounds. The improvement in significance is optimized by a cut | cos^*| < 0.7, 
for which an improvement of 15% is obtained. 
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FIG. 18. Cross sections for — > 77 signal and 77 background as functions of M^^. The 
mass resolution is 0.66-1.0 GeV. 
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FIG. 19. Cross section for 7-quark jet background as a function of M.y«y/, showing the 
rejection obtained by the cuts described in the text. 

Figure 18 shows the Higgs production cross section {(Th) and the cross sections after 
cuts 1-4 and 1-5. It also shows the rate of direct photon background after cuts 1-2, 
1-4, and 1-5. The trigger efficiency for events passing all selection cuts is 98.8%. This 
is included in all signal and background rates. The signals are small compared to the 
background but still statistically significant. 
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FIG. 20: Same as Fig. 19 but for 7-gluon jet background. 




FIG. 21: Jet Background Rates for ^ 77 after various cuts defined in the text. 

QCD jets produce photons from 7r° and r\ decays. These sometimes appear isolated 
and give additional backgrounds to H ^ 77. Both 7-jet events giving 7 "7" and jet-jet 
events producing "77" must be considered. Thus, the jet background to i? ^ 77 is 

au^Y = f^7-jet^(7/jet) + cT2-jet^^(7/jet) , (5.4) 
where 7?.(7/jet) is the probability for a jet to fake a photon. 
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FIG. 22. Drell-Yan e+e^ cross section after event selection cuts (dashed) and lowest order 
direct 77 cross section (solid) . The dots show the rate of fake photon pairs from Drell-Yan 
electron pairs using TZ{'^/e) = 0.005 near the Z and 0.05 elsewhere. 

The probability 7^(7/jet) was determined using a combination of gemf ast for selecting 
isolated events and GEANT for simulating the shower shape cut. QED bremsstrahlung is 
not included, since it is properly included in the higher order QCD cross sections. Since 
7?.(7/jet) is different for quark and gluon jets,^^ separate samples of about 10^ events each 
for and were used to study jet rejection as a function of Mr^ry. Figures 19 and 
20 show the cross sections for the 7-jet background for both quark and gluon jets as a 
function of Mj^^. The shower shape analysis made possible by the fine sampling of the 
GEM calorimeter improves the rejection by about a factor of 4 per photon over isolation 
alone. Since the jet-jet cross section is about 10® times the 77 one, it is not possible to 
generate enough jet-jet events to simulate the background directly. Instead, gg, qg and 
qq events were generated at the parton level, and the probabilities of jets faking photons 
taken from these figures were used to calculate the background. 

Figure 21 shows the cross sections for the jet backgrounds after cut 1 at the parton 
level (including a K-factor of 1.5), after cuts 1-2, 1-3 and after all cuts 1-5. The 7-gluon 
background is reduced to a negligible level compared to the NLL direct photon background. 
The 2 jet background (which is dominated by gluon jets) is reduced to around 15% of the 
NLL direct photon background at 80 GeV. The ^q background, however, remains 40% 
of the NLL direct photon background there. Since a factor of two increase in all the jet 
backgrounds reduces the significance by only 16%, GEM's searches in the ii" — > 77 channel 
are not very sensitive to the QCD jet rejection. 

Drell-Yan e+e" pairs can give a background for — > 77 if the central tracker fails 
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FIG. 23. Higgs mass peaks over subtracted background, obtained with 10 fb \ for i7 — > 
77 searches with Mh = 80, 100, 120, 140 and 160 GeV. Part (a) corresponds to a = 6/8.5% 
in the barrel/endcaps and b = 0.4%, and part (b) to a = 14/17% in the barrel/endcaps and 
6=1%. The GEM photon identification algorithm is used in both cases; see Section 5.1 
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FIG. 24. Higgs mass peaks over subtracted background, obtained with 30fb~^, for H 
77 searches with Mh = 80, 100, 120, 140 and 160 GeV for the GEM detector. 

to identify both charged tracks. Figure 22 shows the invariant mass distribution of Drell- 
Yan e^e^ after event selection cuts and the lowest order direct photon background. The 
Drell-Yan cross section is about 300 times larger than the direct photon background at 



38 




FIG. 25: Number of events and significance for ilf ^ 77 including all acceptances. 

the Z-mass peak, but comparable elsewhere. A GEANT simulation was performed for the 
central tracker to estimate the rate of electrons faking photons, TZ{'-^/e). If one simply 
looks for hits in the central tracker forming a line, the photon acceptance is found to be 
83% for 7^(7/e) = 0.005 and 96% for TZ{'y/e) — 0.05. If one restricts the track to lie in 
a road defined by the pointing resolution of the calorimeter, then the photon acceptance 
improves to 86% for TZ{'j/e) = 0.002. If one uses the vertex found to further restrict the 
road, then a photon acceptance of 91% can be achieved with the same rejection. In what 
follows, we use only the first method, since it allows adequate rejection at all masses, as 
indicated in Fig. 22. 

Figure 23 shows the reconstructed ilf — > 77 mass peaks after subtracting the back- 
ground, for an integrated luminosity of 10fb~\ and Mh = 80, 100, 120, 140 and 160 GeV. 
Figure 24 shows the corresponding peaks after 30fb~^. Figure 25 shows the number of 
signal and background events after all cuts, and the resulting signal significance using a 
±2AMh mass bin. With 10fb~^, GEM could discover a Higgs boson, using this mode 
alone, for 110 ^ Mfj ^ 150 GeV. With 30 fb^^, using this mode alone, it could extend the 
discovery reach down to about 90 GeV and up to 160 GeV. While the heavier masses can 
also be found in the ZZ* mode, the 77 branching ratio is important to distinguish the 
minimal Standard Model Higgs boson from nonminimal ones. 

5.2.2. H(ti/W) 77£X Searches for 80 GeV < Mh < 140 GeV 

The processes Htt — > ^^£X and HW — > ^^£X are complementary to the inclusive 
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H 77, providing essential confirmation. As shown in Fig. 18, the signal cross section 
for the sum of these two channels is of the order of few fb, but the isolated lepton tag and 
photon identification cuts leads to smaller backgrounds, so that the signal-to-background 
ratio is large. In these associated production channels, most of the effects of detector 
resolutions on the reconstructed mass resolution discussed previously remain applicable. 
With a charged lepton in the final state, however, the central tracker is able to determine 
the Higgs vertex unambiguously. 

The main backgrounds to H{tt/W) — > ^^£X searches are: 

1. tt'j^ or 6677 (i^j^j. 

2. qq' — > VF77 ^77, where both photons are radiated from quarks. 

3. qq' ^ (-11 1 where one photon is radiated from the outgoing charged lepton. 

4- QQ/gg Z'j — > £+£"77, where the second photon is radiated from the outgoing charged 
lepton. 

5. qq/gg — > Z'j — > e"7"7, where the fake photon is from misidentification of the electron. 

6. ti — > "77" £X, where both photons are fakes arising from jet fragmentation. 

Since all the higher order QCD corrections to the signal and backgrounds have not been 
computed, leading order cross sections are used for this entire analysis. While tt is initially 
very much larger, all the backgrounds turn out to be comparable after the cuts described 
below. 

Not all of these background processes are included in PYTHIA, so they were calculated 
using a combination of generators. For processes 1 and 2, the initial hard scattering was 
generated using PAPAGENO 3.6."^^ Initial and final state parton radiation, hadronization 
and decays were then generated using PYTHIA. All the final state photon radiation for 
processes 3, 4 and 5 was generated using PYTHIA with the QED radiation option turned 
on. Since the process gg Z'j is not available in either PYTHIA or PAPAGENO, qq — * Z'j 
was increased by 20%.^^ The background (process 6) was generated with PYTHIA with 
QED radiation turned on. The study was carried out with 2.6 x 10^ events. Only 3 events 
survived all cuts; none of them contained two radiated photons. The 6677 background was 
also simulated, and it is small compared to the tt background after isolation cuts. 

The following cuts were made to reject the backgrounds: 

1. \r]i\ < 2.5 and > 20 GeV. 

2. \r]^\ < 2.5 and > 20 GeV. 

3. Photon and lepton isolation with = 5 GeV and R = 0.45 and 0.3 for photon and 
lepton respectively. 

4. Photon and lepton identification by shower shape and track matching. 
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FIG. 26. Signal and background for H{ti/W) —>■ ^^£X after the cuts described in the 
text. 

5. p7p > 40GeV. 

The first two cuts ensured that the lepton and photons could be detected. Cuts 3 and 4 
identify the photons and lepton and reject jet backgrounds. Finally, cut 5 was found to 
reject the backgrounds, especially ti listed above. It reduces the background by a factor of 
2 to 3 while losing only 20% of the signal. The trigger efficiency for events passing these 
selection cuts is 99%. 

Figure 26 shows the Higgs production cross section (an) and the cross sections after 
event selection cuts 1-4 and 1-5 together with the total background rate after all cuts as 
a function of mass. The largest backgrounds are the fake 77 pairs from tt production and 
the QED radiation Z7 £+£"77. A Z-mass veto improves the rejection of the latter 
background, but it does not improve the significance. 

Figure 27 shows the 77 invariant mass spectrum collected in 30fb~^ with Higgs signals 
of 80, 100, 120 and 140 GeV superimposed on the sum of all backgrounds. Figure 28 
shows the significance for the H{ti/W) ^y^lX search with integrated luminosity of 
10 and 30 fb^^. Although the numbers of signal and background events are small, this 
process provides higher significance than the inclusive iiT ^ 77 at very low masses. The 
combination of this channel and the inclusive iif — > 77 provides a more robust search for 
low mass Higgs bosons than either one alone. 

5.2.3. H ZZ* £+£-£+£- Searches for 140 GeV < Mh < 2Mz 

The H — > ZZ^*^ — > decay provides the cleanest Higgs signal. Because of the 
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FIG. 27. Two-photon mass peaks for H{tt/W) — » i'yyX (black area) over sum of all 
backgrounds (gray area), for an integrated luminosity of 30 fb^^ and Mjj — 80, 100, 120 
and 140 GeV. 
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FIG. 28. Number of events and significance for H{tt/W) i'yyX for various masses and 
integrated luminosities. Signal and background are integrated over ±2AMjj. 



four isolated leptons in the final state, most of the QCD background can be rejected by an 
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isolation cut. For a Higgs mass larger than 2Mz, both pairs have an invariant mass 
of Mz, so a mass constraint can also be used to reject background. If the Higgs mass is 
less than 2M^, one (or both) Z-bosons are off the mass shell. In either case, good mass 
resolution is important for rejecting backgrounds. 

The background processes considered are: 

1. zzjzz* ^i+ri+r. 

2. QQZ £+£-£+£- + X, where Q^hoit. 

3. it W^hW^h, in which the two W^-bosons decay semileptonically and the 6-jets give 
isolated leptons. 

The ZZ/ZZ* — > £~^i~i~^i~ background is irreducible. Since the cross section of gg — > ZZ* 
is not yet available in either PYTHIA or IS A JET, its contribution was accounted for by 
multiplying the contribution of qq — > ZZ* by 1.65.^^ 

The following cuts were used to reject the other backgrounds: 

1. \ri^\ < 2.5 and pfp > lOGeV. For electrons, the region 1.01 < < 1.16 was excluded. 

2. Lepton isolation with R — 0.35 and Ej^^ — 5 GeV; see Section 5.1. 

3. Lepton identification and track matching. 

4. 10 GeV < Mjp < 100 GeV and 70 GeV < M^'^^ < 100 GeV to suppress the contin- 
uum background, where M^p and M^j^^ are the low and high invariant masses of two 
£'^£~ pairs. 

The isolation, p^, and mass cuts help reject the heavy flavor background. The trigger 
efficiency for events passing these selection cuts is higher than 98% for the four-electron 
mode and 99% for the other two modes. 

Table 5 shows the cross section after event selections (uaccep) and the corresponding 
background rate for the various H — > £^£~£'^£~ channels. The integrated luminosity is 
10fb~^ While the mass resolution is excellent in the four-electron channel, it is not as 
good in the muon channels because of multiple scattering of the relatively low-energy 
muons in the second muon chamber superlayer. (This is an instance in which simulations 
directly influenced the detector. Changes were made to the muon system design following 
the TDR work to increase its acceptance and resolution for this four- muon process.) The 
signiflcance for all channels combined is also listed in Table 5. Using all the four- lepton 
channels, GEM can discover a Higgs boson in this mass range with an integrated luminosity 
of 10fb~^ The most difficult mass is 170 GeV where the signal is only 5.7cr. Figure 29 
shows the e+e^e+e", ji^ ii~ ji^ ii~ , and £'^£~£'^£~ invariant mass spectra collected in 30 fb~"^ 
for Mff = 140, 150, 160, 170 and 180 GeV superimposed on the sum of all backgrounds. 
The increased integrated luminosity makes even the peak at 170 GeV unambiguous. 
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Table 5. Signal and Background for H ZZ* ^+ /j,- /i- and H ZZ* £+£-£+£- 
for an integrated luminosity of 10fb~^. The signal cross section after cuts is CTaccep- 
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FIG. 29. Higgs mass peaks for H ZZ* t^l~t^l~ over the sum of all backgrounds 
plotted for an integrated luminosity of 30 fb~^, and Mu = 140, 150, 160 and 170 GeV. (a) 
Four electrons, (b) Four muons. (c) All leptons. 



5.2.4. H ^ ZZ ^ £+£-£+£- Searches for 2Mz < Mh < 800 GeV 
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Table 6. Signal and background for H — > ZZ — > for an integrated luminosity 

of 10 fb-^ 
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4.3 


5.3 
600-1200 
1.5 


Background (fb) 


ZZ 


3.0 2.3 


1.0 


.6 


Significance 


Significance 


38 28 


9.7 


4.7 




o 10 




400 



800 400 
M„„(GeV) 



800 



1200 



FIG. 30. Higgs mass peaks and background for H ZZ — > for an integrated 

luminosity of 10 fb"^ (a) Mh = 200 GeV, (b) 400 GeV. (c) 600 GeV. (d) 800 GeV. 

For a heavy Higgs decaying into £'^£~£'^£~, PYTHIA was used to generate the signal 
and the ZZ background. All other backgrounds are negligible. The event selection cuts 
were taken to be: 



1. \t\ < 2.5 and p^j, > 10 GeV. 

2. Lepton isolation with R = 0.3 and E^^ = 5 GeV. 

3. Lepton identification and track matching. 



M 



ZZ 



4M|. 



4. At least one Z with Pt > 

5. \Mu - Mz\ < 10 GeV for both lepton pairs. 
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Cuts 1, 2, 3, and 5 are self-explanatory. If the transverse momentum of the ZZ system is 
neglected, then cut 4 requires sin^* > 0.5 and so reduces the more peripheral qq — > ZZ 
background. The trigger efficiency for events passing these selection cuts is higher than 
99%. 

Table 6 shows the production cross section {(Jh) and the cross section after event se- 
lection cuts ((Taccep) for Mr = 200, 400, 600 and 800 GeV. (A Higgs mass of 800 GeV 
is somewhat greater than the "triviality" bound, but we consider it because it is the 
greatest mass at which the Higgs still appears to give a distinct resonance. Such a mass is 
possible, but it will be accompanied by unspecified new physics at or below a few TeV.) 
The rate of the ZZ background and the combined significance are also listed for an in- 
tegrated luminosity of 10fb~^. For such a data sample, discovery of heavy Higgs bosons 
in the l'^£~l'^l~ modes will be possible for all masses below about 700 GeV. For the 
heaviest masses one would seek confirmation of these four-lepton signals in other channels, 
as we discuss next. We also expect that GEM's capabilities at ultrahigh luminosity (see 
Section 5.6) will permit the discovery of a very heavy Higgs in the four-lepton channel in 
about one year with a data sample of order 100 fb~^. 

Figure 30 shows the £^£^£'^i~ invariant mass spectra collected in 10 fb^^ for Mh = 200, 
400, 600 and 800 GeV superimposed on the background. While the 800 GeV peak looks 
ragged, the calculated statistical significance is 4.7a in 10 fb^^. There is also a systematic 
uncertainty caused by lack of knowledge of the shape of the background. This is discussed 
at the end of the following subsection. 

5.2.5. H ^ ZZ ^ £+£-uu Searches for Mr = 800 GeV 

The branching ratio for H — > ZZ — > i'^i~i'i> is six times that for H — > ZZ — > 
t^t't^l~. Since one Z decays to z^P, the signature is a high-pr Z — > t^t~ plus missing 
energy, and the mass of the Higgs boson cannot be directly reconstructed. The total signal 
cross section is 32 fb. 

The following backgrounds were considered: 

1. qq/gg — > ZZ; the cross section is 73 fb for > 150 GeV. 

2. qq — > Zg and qg — > Zg, with Z — > the cross section is 66 pb for pj, > 150 GeV 
before any cut. 

3. qq, gg — > tt, with t biu; the cross section is 380 pb before cuts. 

To separate the signal from these backgrounds, events with an e+e" or jJi^ jJT pair were 
selected as follows: 

1. |?7^| < 2.5 and > 20 GeV for each lepton. 

2. Lepton isolation with R = 0.3 and E^^ = 5 GeV. 

3. Lepton identification and track matching. 
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Table 7: Signal and Background (fb) for H i+l'i/P 







e^e uu 


Channel 










Signal 


ZZ 


Z + jets 


ti 






16.0 


36.7 


3.3x10^ 


2.0x10^ 


After cuts 1- 


-4 


9.5 


13.8 


1.2x10^ 


53.5 


After cuts 1- 


-5 


7.4 


4.2 


3.1x103 


0.36 


After cuts 1- 


-6 


6.2 


3.2 


2.2 


0.0 






jJL^ jjTvV 


Channel 








Signal 


ZZ 


Z + jets 


tt 


After cuts 1- 


-4 


6.5 


9.4 


8.1x10^ 


36.2 


After cuts 1- 


-5 


5.0 


2.9 


2.1x10^ 


0.24 


After cuts 1- 


-6 


4.2 


2.2 


1.5 


0.0 



4. \M^^-Mz\< lOGeV. 

5. > 250 GeV. 

6. $T > 250 GeV. 

Cuts 5 and 6 are primarily intended to reject the large Z + jets background; see Fig. 31. 
The trigger efficiency for events passing all selection cuts is higher than 99%. Table 7 
lists the cross sections after cuts 1-4, 1-5, and 1-6, for the signal and background in the 
e^e~uu and n^n^uu channels. 

Figure 32 shows the reconstructed transverse mass, M| = 2E^^T (1 - cos A(^), where 
is the transverse energy of the Z and A(j) is the azimuthal angle between the direction 
of the Z-boson and the ^y-vector. This distribution is not sensitive to the degradation of 
the resolution that results from adding to the calorimeter response a 1% nongaussian 
tail with twice the normal width (see Section 4.3). There are 105 signal events over 
a total background of 91 for an integrated luminosity of 10fb~^. Since the signal and 
background distributions are similarly shaped and the signal to background ratio is only 
1.15:1, this cannot be regarded as a convincing discovery channel. A knowledge of the ZZ 
continuum background to ~ 25% would serve to give a 5a systematic-limited significance. 
This knowledge of the ZZ background should be achievable by comparison with WZ and 
Z + jets production. The statistical 4.7o" significance of the four-lepton signal is reduced 
to 4.2cr by the same background uncertainty. Combining the ££££ and ££i7i/ channels with 
Gaussian statistics gives a 6.6a signal for lOfb""^. 

5.2.6. H ^ ZZ ^ (.^(.~33 Searches for Mr = 800 GeV 

The branching ratio for H ZZ £'^£~jj for £ = e and /i is approximately 20 times 
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FIG. 31. Distributions for an 800 GeV H — > £~^£~i'i' signal (solid), ZZ background (dots), 
and Z + jets background (dash) for 10 fb""*^. The background is small, (a) Transverse 
energy of the reconstructed Z. (b) Missing transverse energy {^t)- 

higher than that into all the i~^i~^£'^£^ modes. For Mjj — 800 GeV, the cross section is 
110 fb. The signal-to-background ratio, however, is much worse than in the four-lepton 
channel because of the large background from Z + jets production. In addition, since the 
width of an 800 GeV Higgs is 270 GeV, the signal will be seen only as a broad excess of 
ZZ pairs over the background. 

The following backgrounds were studied in this analysis: 

1. ZW or ZZ with one Z e+e" or ji^ jjT . 

2. Z + jets with Z — > e'^e~ or ji^jjr . 

3. tt with t bW and W ^ ev or jjLv. 

The irreducible ZZ background has a production cross section of 130 fb, for p| > 200 GeV. 
The Z + jets background has a much larger rate of 27 pb for > 200 GeV. The Z- 
mass constraint and high of leptons and jets were used to reduce this background. 
(The possibility of tagging forward jets and vetoing on central jets to enhance the signal- 
to-background ratio is discussed below.) The cross section for the tt background with 
semileptonic decays to e^e~ -\- X and iJi^iJi~ -|- ^ is 380 pb before the Z-mass constraint is 
imposed. A sample of 1.5 x 10^ events was generated for this background. It was largely 
eliminated by pt and isolation cuts and the Z-mass constraint. 

Leptons were selected as follows: 



48 



o 




400 800 
(GeV) 



1200 



TIP-04054 



FIG. 32. Transverse mass Mt of an 800 GeV Higgs signal for H ^ ZZ ^ £+£-iyiy (black 
area) and sum of all backgrounds (gray area), obtained for an integrated luminosity of 10 
fb-i. 

1. jr/^l < 2.5 and pj, > 70 GeV for each lepton. 

2. Lepton isolation with R = 0.3 and ~ 5 GeV. 

3. Lepton identification and track matching. 

4. \Mii-Mz\< 10 GeV. 

5. p^^ > 230 GeV. 

The hadronically-dccaying Z-boson was reconstructed with the following algorithm: 

1. Find all jets at \r]\ < 2.5 using a large clustering cone, R — 0.9. 

2. Find all jets with a small cone, R = 0.3, and match these narrow jets with those found 
using the larger cone. 

3. Require the highest pr jet found with R = 0.9 to have pt > 250 GeV and to be 
composed of two jets found with R = 0.3, each having p^ > 80 GeV. 

4. Reconstruct the mass Mjj of the highest-py dijet with R = 0.9, and require that 



|M. 



J3 



Mz\ < 15 GeV. 



The trigger efficiency for events passing all these selection cuts is higher than 99%. 

Figure 33 shows the reconstructed dijet mass Mjj for the signal and backgrounds, and 
the dijet mass cut is indicated. Note that the background has been multiplied by 0.2 in 
this plot. The mass resolution is about 9 GeV. Table 8 shows the cross sections for signal 
and backgrounds after the lepton and jet cuts for 600 < M^jj < 1000 GeV. 
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Two-jet invariant mass (GeV) 



FIG. 33. Dijet invariant mass distribution for H ^ ZZ ^ for Mh = 800 GeV 

(solid) and sum of all backgrounds (dashed), obtained for an integrated luminosity of 
10fb~^. The background is multiplied by a factor of 0.2. to facilitate comparison of the 
shapes. 

Table 8. Signal and background cross sections (fb) for H — > ZZ — > £'^£~jj for 10fb~^. 
The invariant mass Mnjj signal region is 600 to 1000 GeV. 



Cut 


Higgs 


tt 


Z+iets 


zw/zz 


After lepton cuts 


41 


130 


3240 


20 


After jet cuts 


12 


0.5 


65 


1.5 


Within 600 < M^jj < 1000 


11 


0.3 


42 


0.8 



The invariant mass M££jj of the Z-jet system is shown in Fig. 34. The signal and 
background events passing all these cuts are 110 and 430 respectively for 10fb~^ The 
signal and background have a very similar shape, with no recognizable peak in the mass 
spectrum. Since the signal is 25% of the background, the latter must be known to better 
than 5% to achieve a 5a significance. Thus, this mode cannot be regarded as a discovery 
channel for the Higgs. 

The possibility of tagging forward jets to enhance the signal-to-background was stud- 
63,51 Pqj, this, a segmentation of the forward calorimeter of Arj x A(j) = 0.3 x 0.3 was 
assumed; similar results were obtained with 0.2 x 0.2. A resolution of 100% /\/E © 5% 
was used to smear the hadronic energies. One forward jet was required in 2.5 < |r/| < 4.5 
with E > 2TeV and pT > 50 GeV. After the jet jet cuts described above, the signal was 
reduced to 24 events on a background of 32. Again, the shape of signal and background 
were similar, so there is no significant advantage in tagging a single forward jet. Requir- 
ing two tagged forward jets is expected to reduce the signal and backgrounds to a level 
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FIG. 34. Invariant mass distribution M^+i-jj for 800 GeV H ^ ZZ ^ signal 
(solid) and sum of all backgrounds (dashed), obtained for an integrated luminosity of 
lOfb-^ 

comparable with the four charged lepton channel and, so, will not be of help. Finally, a 
central jet veto was found not to be useful after the high-py cuts on the jets from Z decay 
were imposed. 



5.3. Heavy Flavor Physics 

As we have emphasized, while the characteristic energy of electroweak symmetry break- 
ing lies within reach of multi-TeV hadron colliders such as the SSC and LHC, the scales of 
flavor symmetry and the interactions that break it are unknown; they may be near 1 TeV 
or much higher. Nevertheless, searching for the origin of flavor symmetry and its break- 
down at hadron colliders is just as important as it is for electroweak symmetry breaking. 
Given how little is known of this physics, it is essential that detectors be designed to allow 
for the broadest possible investigations. This can be done only by testing designs against 
a wide variety of scenarios for flavor physics. 

Top-quark physics, in particular, will be an important part of a supercollider's physics 
program for many reasons. Even if the top quark has been found at the Tevatron Collider,^ 
precise measurements of its mass and and other properties may not be possible there 
because of limitations in the detectors and the machine itself. The top quark's apparently 
very large mass of 175 GeV is our most dramatic example of flavor symmetry breaking. 
Thus, if other manifestations of flavor physics are accessible at supercollider energies, the 
top quark will be an important key to flnding them. If charged scalar bosons, H"^, exist, 
then either t — > H^b or — > tb is expected to be an important decay mode. Also, many 
technicolor models contain a color-octet, spin-zero "technieta" boson, tjt, with mass in the 
range 300 — 500 GeV. The ijt would be produced via gluon fusion and is expected to decay 
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into tt pairs. ^^'^'^'^^ Its cross section at the LHC may be enormous. At the same time, 
i-quark production processes can be serious backgrounds to new physics. An example of 
this was encountered in Section 5.2, where it was seen that isolated-photon backgrounds 
arising in tt production were an important background to the ttH'^ £^77 signal for 
production of an intermediate mass Higgs boson. Thus, it is essential to know as much as 
possible about the rate and other characteristics of t-quark production. Finally, the fact 
that t-production gives a pure sample of VF-bosons, with M^^ and mt accurately known, 
may be useful for calibrating the calorimeters for jet energy measurements. 

The work reported in this section was done over a year before the announcement by 
the CDF Collaboration of its evidence for the top-quark. The top-quark masses considered 
in this section are 140 GeV (as generally used in this report) and approximately 250 GeV. 
These bracket fairly well the 175 GeV value reported by CDF. Therefore, we have not felt 
it necessary to repeat the analyses for the CDF central value. This section discusses the 
following examples of top physics in the context of the GEM detector: 

1. Discovery and measurement of the mass of heavy top quarks {mt ^ 250 GeV) pair- 
produced by gg, qq tt and decaying via the standard mode t W^h. In this case, 
one t-quark gives an isolated electron or muon from VF-decay plus a non-isolated muon 
from the decay h — > c/xI7p, and the top mass is determined from the invariant mass 
distribution, -M^^, of the isolated lepton and the non-isolated muon.^^'^^ This study 
may be viewed as a model of using multi-lepton modes to search for a very heavy quark 
decaying into a light one plus a VF-boson. 

2. Discovery and mass measurement of top quarks in tt production, both with mt = 
140 GeV and with mt = 250 GeV. Again, the standard t — > W'^b decay mode is 
assumed. Events are tagged by one isolated high-pr lepton (e or /i) from decay of one 
of the VF-bosons and several high-py jets. The other VF-boson and its parent t-quark 
are observed in their hadronic decay modes. 

3. Discovery and study of a charged scalar H'^ produced in the decay of a heavy top 
quark, t — > H~^b. The scalar will be assumed to decay as t^Pj- one or three 
prongs.^® For this search, one t-quark is tagged using the standard mode t —>■ W~^b — > 
isolated lepton. 

5.3.1. Heavy Top-Quark Detection and Mass Measurement via the Dis- 
tribution 

This analysis^^ uses tt events with two leptonic top decays and an additional non- 
isolated muon from b — > /iX. The top mass is determined from the mass M^^ of the non- 
isolated muon and the isolated lepton of the opposite sign. While the M^^ distribution is 
broad, this analysis avoids hadronic energy measurements and minimizes the sensitivity to 
the top production dynamics. 

IS A JET 6.36 with the EHLQ-1 distribution functions was used to generate 60K tt 
events for each of the masses mt = 200, 230, 250, 270 and 300 GeV. (The t-quarks 
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were generated with 50 GeV < pt < 1000 GeV, and the appropriate decays were forced.) 
The signal events are those with one isolated electron and one isolated muon of opposite 
sign, together with one non-isolated muon of either sign. Choosing isolated leptons of 
different flavors eliminates backgrounds from — e~'"e~, n^fi~. The signal cross sections 
(assuming B{W iu) = B{b ^ + X) = 1/9) for mt = 200, 250 and 300 GeV were 
found to be 19.9 pb, 8.2 pb and 4.0 pb, respectively. 

The signal cross sections are much larger than all backgrounds, as discussed below. 
To optimize the separation between isolated leptons from VF-decay and non-isolated ones 
from 6-decay, the following cuts were imposed to select isolated electrons: 

\Ve\ < 2.4 
J2 Et{EC) > 40 GeV 

(5.5) 

J2 Et- Et{EC) < 10 GeV 

R=OA R=0.2 

0.8 < Et{R = 0.2)/pr(CT) < 1.2 , 

where Et is the transverse energy in a cell of the full calorimeter; i?r(EC) is the transverse 
energy in EM calorimeter cells; pr(CT) is the electron momentum as measured in the 
tracker. As usual, the sums are over rj — (f) cones of radius R. The cuts for isolated muons 
were: 

\V^\ < 2.4 

PT>40GeV ^gg^ 
J2 Et- Y Et< 2GeV. 

Non-isolated muons were required to satisfy: 

I^mI < 2-4 
PT>20GcV 

Y Et- Y Et>IQ GeV . 

R=QA R=Q.l 



The choice R = 0.2 for the cone defining the electron energy is not critical because 
the electron energy measurement does not need to be very precise here. The and 
M^^ distributions turn out to be very similar for this choice of isolation parameters and 
not very sensitive to moderate changes. For the relatively low-py muons in the top-quark 
decays, the muon energy may be corrected simply by adding the average energy loss in the 
calorimeter for a given momentum and rapidity. For a top mass of 250 GeV, the isolation 
cut (5.6) retained 78% of the muons from VF-decay while rejecting 99% of the muons from 
6-decay. The non-isolation cut in (5.7) accepted 89% oi h ^ ji + X and rejected 91% of 
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FIG. 35. Mifj^ distribution for an isolated lepton and a nonisolated muon of opposite sign 
for mt = 200 (solid), 250 (dotted), and 300 GeV (dashed). 

There was a 30% contamination of the ti — > e^^^ /xj^^^j + /Xnon-isol signal due to non- 
isolated muons in the decay of c-quarks which, in turn, came from the decay of the "wrong" 
t-quark or from g cc. This background was reduced to the 15% level by selecting events 
in which the isolated lepton and non-isolated muon of opposite sign are close to each other 
and, hence, more likely to be from the decay of the same t-quark. Events were required to 
satisfy 

5(l){£fj,) < 90° . 

The acceptance of this cut was found to be 67%, 56% and 49% for mt = 200, 250 and 
300 GeV. A detailed simulation showed that the reconstruction efficiency was greater than 
95% even for non-isolated muons from t ^ b ^ decays. The overall acceptance ranged 
from 1.1% for mt = 200 GeV to 1.8% for mt = 300 GeV both for the tj^ ej^i /"non-isol 
events and for the //J^j iJ,^^^_ 

-isol 6V6nts. These efficiencies are low because the analysis was 
designed to obtain a very clean sample. The total number of £^^^ A*non-isol ^'^^^^^ expected 



with an integrated luminosity of lOfb 
300 GeV. 



-1 



is 4500, 2400 and 1400 for mt = 200, 250 and 



Backgrounds to the ti — > e^^^j /x^^j Afnon-isol signal from production of + jets and 
+ jets were considered. The most important contributions come from 6-quark jets 
and — > T+r" + fj/^ + ^t- After the selections described above, the + jets 

background was found to be less than 9% of the signal for the 250 GeV top quark; the 
Z^ + jets background was negligible. 

The shapes of the Mg^ and M^^ distributions are nearly identical. Thus, adding them 
together to form a Mf^ distribution adds statistical weight to the mt determination without 
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FIG. 36: Mean value of M^^ vs. m^. 

introducing significant systematic error. The M^^j distributions arc plotted in Fig. 35 for 
several values of nit- The mean value of M^^ should vary linearly with nit and does, as 
can be seen from Fig. 36(also see Ref. 69). A statistical error of 2.3 GeV on the mass of a 
250 GeV top quark is expected for 10 fb^^ of data. 

The systematic error in this method of determining mt arises from imperfect knowledge 
of the 6-quark fragmentation function and of the px distribution of the t-quarks. These 
affect the momentum distributions of the non-isolated muon and the isolated lepton, re- 
spectively. The effect of heavy-quark fragmentation was studied by varying the parameter 
e in the Peterson fragmentation function from its nominal value of e = 0.006.^ The range 
considered, la = 0.002, led to a variation of 2.3 GeV in the mass of a 250 GeV t-quark. A 
measure of the systematic error due to the pj'(^)-distribution was obtained by varying the 
amount of initial state radiation in tt production. This resulted in a 3.4 GeV change in the 
mass determined for the 250 GeV t-quark. Finally, next-to-leading-order QCD corrections 
to the tt cross section increase the magnitude of the cross section by about 60%, but do not 
significantly change the shape. ™ Still higher-order corrections are expected to amount to 
±15%. We conclude that, with an integrated luminosity of 10 fb^^, it would be possible 
to use the M^^ method to determine the mass of a 250 GeV top quark to within an error 
of ±2.3 GeV (statistical) and ±4.1 GeV (systematic). 

5.3.2. Top-Quark Detection and Mass Measurement via the Mjjj Distribution 

The most direct measurement of the top-quark mass comes from the nonleptonic decay 
t — > Wb — > qq'b, giving three jets. This is also the most precise measurement if systematic 
uncertainties associated with jet definition and energy measurement are under control. 
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Furthermore, it is important for flavor physics spectroscopy to observe the top quark in 
nonleptonic decay modes. 

ISAJET was used to generate 300K tt events for mt = 140 and for 250 GeV.^'^ As above, 
the top quarks were generated with 50 GeV < pt < 1000 GeV. Events were selected in 
which one W decayed to an electron or muon while the other W decayed nonleptonically. 
The physics background to tt i^^^ + jets comes mainly from production of W + jets. 
Since both the signal cross section and the signal/background ratio are much higher than 
at the Tevatron,^ it is possible to make cuts that render this background unimportant. 
Thus, a non-isolated muon tag was not required. The ISAJET cross sections for these 
events, including VF-boson branching ratios, are 4.1nb and 0.44 nb for the two values of 
mt- The isolated leptons were required to satisfy the cuts in Eqs.(5.5) and (5.6). 

To eliminate combinatorial backgrounds in the multijet mass distributions, the t-quark 
and individual jets were forced to be at high pT and to be in the hemisphere opposite the 
isolated lepton. For such events, the three jets from the decays of the t tend to be close, and 
a cone of small radius R = 0.4 was therefore used. Individual jets were required to have 
Pt > 30 GeV (50 GeV) for mt = 140 GeV (250 GeV). In both cases, jets were required 
to satisfy S4>{£jet) > 90°. Finally, the three highest-pj- jets satisfying these constraints 
were required to have |pT(3jets)| > 200 GeV (300 GeV) for mt = 140 GeV (250 GeV). The 
combined geometrical acceptance and efficiency of these cuts was found to be 0.70% and 
1.0% for mt = 140 and 250 GeV. The number of events obtained per 10 fb^^, taking into 
account trigger and reconstruction efficiencies and these acceptances, were 240K and 40K 
respectively. The W + jets background was found to be negligibly small for the 140 GeV 
case. It contributed 15% to the three-jet mass spectrum for 250 GeV, but this had no 
important effect on the determination of mt- 

Consider first the case of mt = 140 GeV. The W — > 2 jets and t — > Wb — > 3 jets 
mass distributions, Mjj and Mjjj, are shown in Figs. 37 and 38. Here, the trijet search 
region was 110 GeV < Mjjj < 170 GeV, and at least one dijet pair was required to satisfy 
\Mjj — Miy\ < 20 GeV. All dijet pairs passing the cuts appear in Fig. 37 (a). The 
fitted VF-mass peak is at 75.8 GeV in this figure, and the resolution is 13.9 GeV. The 
large combinatorial backgroimd is a consequence of the kinematics: for mt = 140 GcV 
and My^/ = 80 GeV, the three jets tend to be roughly equidistant in 77 — ^ space. The 
background comes from picking up the 6-jet from the same or, less frequently, the other 
t-quark decay. The signal-to-noise ratio in the W^-peak region is about 2:1. 

This signal-to-noise can be improved considerably by requiring that one of the jets 
passing the above cuts be tagged as a 6 jet by the tracker. Figure 37 (b) shows the Mjj 
distribution for events with an identified b jet which is then excluded from the dijet mass. 
The b-tag efficiency for this plot was assumed to be 20%. The W signal-to-noise ratio 
is improved to 3.9:1. The fitted W mass in this plot is at 76.9 GeV with a resolution of 
8.4 GeV. We emphasize that this b tag is not required for determining mt- This is seen in 
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FIG. 37. Dijet mass distribution in tt events with mt = 140 GeV. (a) No 6-tagging was 
assumed, (b) 6-tagging was assumed with 20% efficiency. 



Fig. 38, where the top-quark peak appears clearly above little combinatorial background. 
The mass was determined to be 138.2 GeV with a resolution of 8.1 GeV. 

For nit = 250 GeV, the dijet and trijet mass distributions are shown in Figs. 39 and 
40. The search region was 200 GeV < Mjjj < 300 GeV. As above, at least one dijet 
was required to have invariant mass within 20 GeV of M\y and the invariant mass of the 
two closest jets was plotted. There was no need for a 6-tag to sharpen the dijet mass 
distribution, since for such a heavy top-quark the two jets from the W are closer to each 
other than either is to the b jet. Thus, the combinatorial background under the VF-peak 
is much smaller if one selects the closest two jets. The fitted W mass was found to 
be 79.8 GeV, with a resolution of 7.1 GeV. The top-quark mass was determined to be 
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FIG. 38: Trijet mass distribution in tt events with nit = 140 GeV, without 6-tagging. 



3000 



° 2000 



> 
(5 



> 

LU 




1000 



80 120 
Mjj (GeV) 



200 



TIP-03962 



FIG. 39: Dijet mass distribution in tt events with mt = 250 GeV. 

247.4 GeV with a resolution of 14.7 GeV. 

The statistical uncertainty on the top-quark mass for 10fb~^ of data is approximately 
0.03 GeV (0.11 GeV) for mt = 140 GeV (250 GeV). Systematics dominate the error in this 
mt measurement. The largest effect is the uncertainty in the jet energy measurement. 
More detailed studies of the calorimeter — beam tests as well as simulations — will be 
required to yield adequate correction functions. It will be helpful that the t Wb process 
is self-calibrating: The position and width of the W peak will help calibrate the energies 
of jets in t-decay and determine the systematic uncertainty in mt- 



The top quark will be a signal of flavor physics. One example is the rjT, a spin-zero 
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FIG. 40: Trijet mass distribution in tt events with mt = 250 GeV. 

color-octet boson occurring in many technicolor models. ^^'^^ In modern models, the tjt 
is expected to have a mass in the range 300 - 500 GeV. ^^-^^ It is produced in pp collisions 
via gluon fusion with a cross section of order 1-50 nb, depending on its mass, the number 
of technicolors, and other model-dependent factors. It is expected to decay predominantly 
to ti with a width of order 10-100 GeV, again model-dependent. 

In a search for the rjT, one would tag ti production in the ^j^^jj-H jets mode and look for an 
enhancement in the ti invariant mass distribution. This invariant mass can be determined 
up to a quadratic ambiguity by assuming that the missing- is due to — > £^z/^. The 
resolution on this invariant mass for M^j ~ 400 GeV is about 15 GeV plus the contribution 
from the resolution. For a 20 GeV wide tjt, and assuming that the M^j resolution is 
30 GeV, the ijt appears as an enhancement 75 GeV wide. The underlying ti cross section 
in such a region is about 3 nb for mt = 140 GeV and M^^ = 300-500 GeV. There should 
be no difficulty discovering the tjt over such a background. 

5.3.3. Discovery of a Charged Scalar in the Decay of a Heavy Top Quark 

Charged color-singlet scalar bosons, i?^, occur in multi-Higgs-doublet models, in all 
supersymmetric extensions of the Standard Model and, typically, as technipions in models 
of dynamical electroweak symmetry breaking. Generally, their couplings are Higgs-like. 
That is, tend to decay to the heaviest fermion pairs kinematically accessible and they 
often must be in the same generation. Since are color-singlets, they are copiously 
produced only if there is a heavier quark which can decay into them. 

We study discovery of an in the decay products of a heavy top quark. The masses 
mt = 250 GeV and Mjj+ = 150 GeV were used for the discussion that follows. At the end, 
results for other mass combinations will be described. It was assumed that both t — > W'^b 
and t — > H^b decays are allowed. The charged scalar was assumed to decay into cs and 
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T'^i'r- Only the decay mode r^Ur is considered here,^^ since for H'^ — > cs, it should 

be straightforward to discover the and measure its mass using the same technique as 
described for the study of t ^ W~^b 3 jets in the previous section. 

The various branching ratios for t and decay are model-dependent. In the absence 
of experimental support for any particular model, it is appropriate to assume only that 

B{t W+b) + B{t H+b) = 1 and B{H+ cs) + B{H+ t+v^) = 1, and to study 
the reach of the detector for as a function of the branching ratios. For comparison 
with other simulations J^'^"^ we also present results in terms of the two-Higgs-doublet model 
occurring in the minimal supersymmetric standard model (MSSM).^^ 

The presence of t ^ H^b, followed by T~^Ur, is signalled by a breakdown of 

lepton universality expected if only t W'^b were allowed (also see Refs. 72 and 73). In 
addition to the tt+z^t decay mode used as a r-tag in those references, all hadronic r-decay 
modes were used here. 

ISAJET 6.50 was used to generate 15K each of ti W+W^bb and tt WH^bb 
events. The ISAJET decay table was updated to include all major r decay modes. 
Hadronic modes account for 64.5% of all decays while the tt^ mode is only 12%. Sig- 
nal events were selected by requiring an isolated electron or muon and a tau candidate. 
The criteria for an isolated fi^ were taken to be: 

I^mI < 2.4 
PT > 20 GeV 
J2 Et- ^ Et<2 GeV 

R=QA i?=0.1 

Muons were required to be fully reconstructed according to the efficiency parameterized 
in gemf ast. For an isolated the criteria used were 

\rie\ < 2.4 

> 50 GeV 

i^=o.3 

^ Et{^C) > 40 GeV 

7?=0.2 

^ Et- Et(EC) < 10 GeV 
0.8 < Et{R = 0.2)/j9r(CT) < 2.0 . 

For analysis of the r-lepton's two-body decays, it is important to account properly 
for the polarization it has in and decays. The decay — > r^Vr conserves 
chirality (the same as helicity for a high-energy r), while — > t^u-j- maximally violates 
it. Then, for example, the tt^ occurring in r"*" — > tt^Ut tends to follow the direction 
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(5.8) 



(5.9) 



of motion and so has higher pT than tt"*" from W 



+ 



decays. This effect is enhanced 



here because is larger than M^r. The polarization correlations were implemented in 

ISAJET for two-body r decays. The r-polarization effects are small for the three-prong 
decays and require no special simulation. 

Isolated r candidates were selected by requiring: 

hrl < 2.4 

iVch = 1 or 3 

Et>50 GeV 

R=0.3 

S4){eT) > 100° 

J2 Et- Y Et< lOGeV 



(5.10) 



R=QA 



R=Q.2 



For the charge multiplicity (A^ch) cut, tracks were required to have pj- > 1 GeV and to 
lie within a cone of = 0.1 around the calorimeter jet axis. The leading track in this 
cone was required to have pj- > 15 GeV. Electrons from r decays were not included in 
the r sample if they passed the isolated e criteria in Eq.(5.9). Muons from r-decays are 
predominantly isolated and are thus rejected by the absence of significant activity in the 
calorimeter. The central tracker was not used to tag r leptons by their displaced vertices. 
If it were, it would further enhance the significance of the nonuniversality signal. 

In addition, we required a h jet tagged by the central tracker. The 6-jet criteria were 
|r7b| < 2.4, d(j){lh) < 100°, and scalar Et > 30 GeV. A tagging efficiency e(6-tag) = 20% 
was assumed in this study. Alternatively, b jets could be tagged by a non-isolated muon, 
but this would decrease the signal significances somewhat. 

The number of standard deviations A'^^r by which the number of r events exceeds the 
expectation from lepton universality is 



N{WH^ ^ 4o1t) 



y/N{WW ^ 4oit) + N{WH± ^ 4oit) ■ 
Here, the event numbers are given by 

N{WW 4o1t) = '2Nj^B^{t W+b) 

X B{W Ivi) B{W Ti^r) B{t hadrons) 

X €{W £) e{W t) e(6-tag)e(t-tag) e(t-trigger) ; 

NiWH"^ 4o1t) = 2A|^ B{t^ W^h) B{t ^ H+b) 

X B(W iu^) B(H+ TVr) B{t hadrons) 
X e{W £) e{H+ T)e{b-tag) e(t-tag) e(t-trigger) 



(5.11) 



(5.12) 
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FIG. 41. Significance of the t — > H'^b, — > t'^Vj- signal vs. B{H^ — > T'^Vr) for various 
5(t ^ H+b). Here, = 250 GeV and Mh+ = 150 GeV. 

For a 250 GeV top quark, the number of events produced with an integrated luminosity of 
10 fb^^ is A^^^ = 1.5 X 10^. The efficiencies e{W r) and e{H^ r) are the ratios of the 
numbers of true r-leptons passing the above cuts to the number generated, including r — > e 
events misidentified as hadronic r-decays. For rrit = 250 GeV and = 150 GeV, these 

efficiencies were found to be e{W — > r) = 9.8% and e{H^ — > r) = 14.7%. The efficiency 
for finding W — > ivi events was 36% for electrons and 45% for muons. The average value 
e{W — > = 40% was used. The top-tagging efficiency, e(t-tag) = 69%, is the fraction of 
top events remaining after the tau-selections. The efficiency for triggering on top quarks 
at Level 1 was found to be 93%.^^ 

The backgrounds to the tt W^H"^ £^ + t^ + X signal come from (1) ti WWbb 
and WH'^bb events in which and i?^ decay to jets which fake a r; (2) W + QQ 
events with Q = c,b; and (3) bb production. To study the first background, 40K two-jet 
events were generated with in the range 50 — 1000 GeV. All jets found by gemfast 
with Et > 25 GeV were subjected to the r-selection criteria above. A rejection factor of 
7^(r/jet) = 0.0027 was found. Since there are two to three jets in the hemisphere opposite 
the isolated lepton, this background amounts to 8.4% of the tt — > WW ^igoi -l- r signal 
and 1.4% X (1 - B{H+ tv^))/B{H+ tv^) of the tt WH + r signal. 

Including these backgrounds reduces by 3%. The other backgrounds have been shown 
to reduce AT^ by less than 3%.'^^'^^'^^ 

The significance for 10fb~^ is plotted against B{H^ — > t'^Vt) in Fig. 41 for B{t — > 
H^b) = 0.01 — 0.99. Backgrounds are not included in these plots because they are so 
small. A signal {Na > 5) for the t — > — > r"*" decay channel can be discovered if 
B{t H+b) > 0.02 and B{H+ r+i^r) Z 0.05. 
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FIG. 42. Significance of the t — > H'^b, H'^ — > r'^i'j- signal vs. tan/? for various mt = 
250 GeV and Mjj+ = 225 GeV (dash-dotted); 250 and 150 GeV (dotted); 150 and 125 GeV 
(sohd); 150 and 100 GeV (dashed). 

We have also computed N„ for the two-Higgs-doublet model used in the MSSM. In this 
model, B{t — > H'^b) and B{H^ — > t'^Vt) are determined by a single parameter, tan/?, the 
ratio of the vacuum expectation values of the two Higgs doublets. The significances were 
determined for four cases: {mt, Mh+) = (250, 225), (250, 150), (150, 125), and (150, 100) 
to facilitate comparison with Ref. 72. The significances expected in GEM are shown in Fig. 
42. For all but the heaviest mass combination, GEM is able to detect a 5cr enhancement 
in T production for tan/? > 0.4. For {mt, Mjj+) = (250, 225), the discovery level covers 
the range 0.7 < tan/? < 2.0 and tan/? > 20. By tagging the r lepton in all its hadronic 
modes instead of just r — > ttUt-., the significances have been increased by a factor of 2.3 — 5. 

5.4. Jet Physics at Large Transverse Momentum 

Although GEM emphasized precision measurement of photons and leptons up to the 
highest SSC luminosities, the accurate identification and measurement of jets is also im- 
portant. The search for a high mass Higgs boson (Section 5.2) and heavy fiavor physics 
(Section 5.3) are two examples of new physics requiring accurate understanding and mea- 
surement of jets. Hadronic jets can also be important backgrounds to new physics; the 
search for H'^ 77 is an example of this. The errors in jet measurements introduced 
by analysis effects, such as jet definition and clustering algorithms, and by instrumental 
effects, such as detector resolution, e/h for the calorimeters, and cracks and dead spaces, 
must be carefully studied and kept under control. This section discusses the main issues 
involved in measuring the high-py jet cross section accurately. The search for signs of 
quark substructure in deviations from QCD of the jet cross section forms the context for 
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this discussion. 

If quarks and leptons are composite, with structure at the scale A, the most visible 
manifestation at subprocess energies V§ <^ A is the presence of four-fermion contact 
interactions, C\, involving the composite quarks and leptons. These interactions induce 
terms in the cross section for dijet and dilepton production that are of order tts/A^, leading 
to significant excesses at "low" s. It is known from experiments at e+e^ and hadron 
colliders that A ^ 1 — 2 TeV, above the scale of electroweak symmetry breaking. ^'^^ Thus, 
C\ must be SU{3) ® SU{2) ® U{1) invariant, and the composite quark and lepton fields 
appearing in it are electroweak, not mass, eigenstates. This raises the possibility that 
unacceptably large flavor-changing neutral currents will appear in the contact interactions. 
The most stringent limit on such interactions comes from the allowed magnitude of new 
effects in the neutral kaon system. If |AS'| = 2 contact interactions exist, they must have 
A ^ 500 TeV, well above the reach of the SSC. For the discussion of contact interactions 
accessible at the SSC, therefore, we shall assume that C\ is symmetric under interchanges 
of the three generations of quarks and leptons. 

Quark substructure shows up directly as an excess of jets at high pT or v^. In this 
section the contact interaction 

All 

-^qq = Qhal^Qha QuIij-QlI, (5-13) 

is used as a model to modify QCD jet production. Here, Qlo = i'^a^ daji are left-handed 
quark fields and a,b = 1, 2, 3 label the generations. This model for the four-quark contact 
interaction is essentially the one discussed in Ref. 11 except that all quarks are taken to 
be composite. PYTHIA 5.6 was used to generate the jet events for QCD and the quark 
substructure signal.^* Several different choices of parton distribution functions (PDFs) 
were used and results compared: EHLQ Set 1,^^ the CTEQ Set IL,^^ and Morfin-Tung 
Set 2.^9 The signal region lies above Et = 4 TeV for the jets. The study described here 
was made for an integrated luminosity of 10fb~^. We shall show that GEM could reach a 
substructure scale of about 25 TeV with this data sample. 

The most important issue in searching for quark substructure is to be certain that 
an observed excess of high-i?r jets is not an artifact of the detector nor of the analysis. 
The jet cross section must be normalized to the QCD expectation at lower Ej^ to reduce 
uncertainties due to luminosity (approximately 3% at the Fermilab Tevatron^). Since the 
cross section for jets at the lowest transverse energies is not well-known theoretically, the 
normalization region used is the middle of the jet Et spectrum. Finally, a scheme must 
be developed that corrects for the calorimeter's lack of compensation. 

For this study signal events were generated with total scalar transverse energy, ^ Et ^ 
9 TeV and the normalization sample was generated with J2^T > 4.8 TeV. The only 
kinematic cut imposed on jets is \r]\ < 1.1. This enhances the roughly isotropic signal 
relative to the forward-backward peaked QCD background. More than one jet can be 
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FIG. 43. Inclusive jet spectrum for substructure at the scale A = 25 TeV (dashed) and 
for QCD (solid). 

taken from a single event. For these central high-E'y jets, the triggering efficiency is close 
to 100%. The discovery criterion adopted for this analysis is an excess of 100 events in an 
£'r-region in which the observed cross section is twice as large as the QCD expectation. 

Jet clustering was done with a fixed-cone algorithm using a large cone, R = 0.9. The 
large clustering radius was chosen to reduce energy loss out of the clustering cone. This 
cone algorithm was found to be efficient and insensitive to detector variations and to have 
good angular resolution.^^ If the jet-cone radius was decreased to 0.7, the main effect was 
to shift the energy scale of the jet-^^j- spectrum downward by 1.2%. This shift should be 
calculable in perturbation theory. 

A full experimental analysis would have to include the development of a jet energy 
correction function and an unsmearing procedure for the inclusive jet-E'^ spectrum. The 
correction function must include the effects of the underlying event, showers spreading 
outside the clustering cone, and detector noise, inefficiencies and nonlinear ities. Most of 
these corrections are reasonably well-understood and tend not to be important for very 
high-i?y jets. The most important problem is to determine the jet energy scale. Jets 
whose Et is mismeasured upward easily produce a false compositeness signal. The jet 
energy scale for the calorimeter would be established using test beam and collider data 
and the energy scale correction would depend on, among other parameters, the jet rapidity 
and the fraction of electromagnetic energy observed in a jet. 

The reach in A is limited by incomplete knowledge of the parton distributions and 
of the jet energy resolution and scale. A number of jet energy reconstruction schemes, 
particularly weighting algorithms which attempt to boost the hadronic part of a shower, 
have been used to improve jet energy resolution and linearity.^^'^^ We have estimated 
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FIG. 44. Number of jets expected with 10 fb ^ with Et > E^^, for various parton distri- 
butions. The sohd hue is QCD (no substructure), dotted hne is A = 30TeV and dashed 
line is A = 25TeV. 

the systematic effects due to nongaussian tails in the energy response and to nonlinearity 
arising from lack of compensation in the calorimetry. Nongaussian tails were modeled in 
gemf ast by adding a second Gaussian of 1% the amplitude and three times the width 
of the normal Gaussian energy resolution, as described in Section 4.3. To estimate the 
uncertainty due to parton distributions we compared the rates of the CTEQ-IL, EHLQ-1 
and MT-2 distribution functions. 

The inclusive cross sections for jets with |?7| < 1.1, using the CTEQ-IL distribution 
functions, are shown in Fig. 43 for A = 25 and QCD (A = oo). Fig. 44 shows the number 
of jets expected in 10fb~^ with Ej- > E™^. It is clear that a quark substructure signal 
at the scale A = 25 TeV could be discovered easily in this integrated luminosity with a 
detector whose calorimeter is as linear as the one modeled here. The variation in jet rates 
obtained using the CTEQ-IL, EHLQ-1 and MT-2 distributions is indicated in Fig. 44. 
The differences are small compared to the signal and produce a 10 — 20% variation in the 
observable A according to our discovery criterion. 

Nonlinearity in the charged pion response was modeled using a GEANT simulation of 
the GEM detector incorporating the expected e/h values for each calorimeter. Figure 45 
shows the charged pion response, averaged over \rj\ < 1.1, as a function of energy. The 
relative response has been set to 1.0 at 200 GeV. The solid curve is the average response 
obtained using energy-independent weights for each calorimeter layer; the dashed curve 
shows the response with energy-dependent weighting. The extrapolation to high energies 
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FIG. 45. Calorimeter response to single pions relative to that at 200 GeV, using energy- 
independent weights (solid curve) and energy-dependent weights (dotted curve). 

Table 10. Number of jets expected with 10fb~^ using different energy- weighting schemes 
as described in the text, and the CTEQ-IL parton distributions. 
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shown here gives noticeable errors in both cases, but mimics the extrapolation from test 
beam energies that will have to be carried out when calibrating the real detector. 

The effect of the nonlinearity in the two weighting schemes is illustrated in Table 10. 
The CTEQ-IL distribution functions were used here. As expected from Fig. 45, the energy- 
independent weighting scheme tends to increase the measured Et of the jets, while the 
energy-dependent scheme reduces jet-ET- The effects are small, of order 5% and the scale 
A = 25TeV is still easily within reach for a data sample of 10fb~^. 

Finally, we can estimate the reach in the quark compositeness scale A attainable with a 
data sample of 100 fb~^. We anticipate no special difficulties for high-E'j- jet measurement 
associated with operations at £ ~ 10^^cm~^s~^. Thus, the reach in A can be determined 
from the fact that the subprocess cross section goes as s/A^. This yields A ~ 45TeV, a 
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factor of 40 greater than the hmit set by existing hadron coUider data. 



5.5. Supersymmetry 

Supersymmetry (SUSY) is theoretically attractive because it eliminates the quadratic 
divergences in the Higgs sector and so allows light elementary Higgs bosons to occur 
naturally. Its study also provides a good testing ground for many aspects of GEM, including 
missing energy, jets and leptons. The minimal supersymmetric standard model^^ has two 
Higgs doublets and superpartners (denoted by a tilde) for all normal particles. In particular 
there are four neutralinos, x^, which are linear combinations of the partners of the photon, 
Z, and neutral Higgs bosons, and two pairs of charginos xf- There is a conserved R parity 
carried by all superparticles, so they must always be produced in pairs and decay to the 
lightest supersymmetric particle, which is absolutely stable. We shall assume that this 
lightest particle is the lightest neutralino Xi- If SUSY is broken at the electroweak scale, 
the masses of all of these particles should be less than about 1 TeV. 

In this section we shall discuss only searches for gluino and squark production. Neu- 
tralino, chargino, and slepton production are all important parts of the search for su- 
persymmetry but have smaller cross sections. While MSSM Higgs bosons have rather 
different production cross sections and decay modes than Standard Model ones, the meth- 
ods of searching for them are generally similar. Figure 46 shows the regions of MSSM 
parameter space in which GEM would be able to detect one or more of the neutral Higgs 
bosons using the 77 and £'^£~£~^£~ decay channels. The r+r" decay modes are not shown 
but appear to offer potential. 

5.5.1. Signature for Gluinos and Squarks 

Since the lightest supersymmetric particle Xi is neutral and interacts weakly with 
matter, it escapes from the detector. Thus, one of the basic signatures for SUSY is missing 
transverse energy, plus multiple jets. A stringent test for GEM's missing energy 
resolution is to be able to detect, in this mode, gluinos and squarks with masses as light 
as 300 GeV. This is near the limit expected from the Tevatron and is also the mass range 
expected in some SUSY grand unified models. ^^'^^ The MSSM typically produces cascade 
decays from one supersymmetric particle to another. The events can have many jets and 
leptons, and the missing energy from the final lightest supersymmetric particle Xi can be 
small compared to the parent mass. A typical decay sequence for a relatively light squark 
and gluino with M~ > M~ might be 

g^xtud, (5.14) 

— h ~Q + 

Decay chains can be even more complex for heavier masses. All of these possible decays 
are included in the version of ISAJET^^ used for this analysis. 
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FIG. 46. Regions of the {Ma, tan/3) plane in which GEM could discover the neutral Higgs 
bosons of the MSSM at the 5a significance level, in the 77 and decay channels, 

with integrated luminosities of 30fb~^ and 100 fb^^. Also shown is the region covered by 
LEP II at ^/s = 190 GeV with an integrated luminosity of 500 fb~^. 

Table 9. Choices of MSSM parameters for the three cases considered. These were chosen 
to have different event topologies and to span the whole mass range. All the squarks and 
all the sleptons are taken to be degenerate for simplicity. All masses are in GeV. See 
Ref. 52 for the notation. 



Parameter Case I Case II Case III 



M~ 300 350 2000 

M~ 600 325 2500 

Mj 500 200 1500 

Ma 300 300 300 

-300 -300 -1000 

tan/3 2 2 2 



There are a number of other parameters in the MSSM besides the gluino and squark 
masses, and it is beyond the scope of this study to explore the MSSM parameter space 
completely. Instead, the representative choices listed in Table 9 have been considered. 
Case I has a light gluino and a heavier squark; it is generally similar to the models of 
Ref. 81 and to the case considered in previous GEM studies. Case II has a squark 
slightly lighter than the gluino and is generally similar to the models of Ref. 82. Since 
g ^ qq dominates for M~ > M~, the signatures in this case are similar to those for squark 
pair production. One might think that this case would be more difficult to detect because 
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the events contain just two hard jets from q — > x^qf. It is actuaUy easier, because the 
branching ratios for qi — > x^"? cind 5l turn out to be large and to provide multijet 

signatures, and the dominant decay — > Xi? gives a harder distribution. Finally, 
Case III has all the masses pushed close to the upper limit for SUSY to be related to the 
electroweak scale. It tests the ability of GEM to cover the top of the plausible mass range 
for weak-scale supersymmetry. 

For the three cases, samples of 70K, 25K, and 35K, respectively, of gluino and squark 
signal events were generated with a version of ISAJET containing all the MSSM decay 
modes. The total production cross sections for all combinations of gluinos and squarks 
are 

Case I: (7 = 8.27nb, 

Case II: (7 = 7.60nb, (5.15) 
Case III: a = 0.81 pb . 

The Monte Carlo statistics are therefore small compared to those obtained in 10fb~^ for 
the first two cases but comparable in the third. This is reflected in the error bars on the 
plots shown below. 

The signal events are characterized by multiple jets and large missing energy. For 
the lower masses in Cases I and II the dominant Standard Model physics background 
comes from heavy flavors decaying into neutrinos, and the dominant detector-induced 
background comes from mismeasuring QCD jets. A total of 1.5M QCD jets of all types 
in ten px ranges covering 50 < pt < 3200 GeV was generated with ISAJET to determine 
both kinds of backgrounds. For the high masses in Case III, the backgrounds from W ^ lu 
and Z ^ vv are also significant. A total of 40K W ^ iv and 80K Z ^ vv events was 
generated covering the same pt range. 

The detector response to all events was simulated with gemf ast. The missing energy 
was calculated using the single-particle pj- resolution of the forward calorimeter determined 
from GEANT plus an additional 1% nongaussian tail three times as wide as the main peak, 
as described in Section 4.3. The effect of this nongaussian tail is small compared to the 
effects of angular resolution in the forward calorimeter and of the hole for the beam pipe, 
so its exact parameterization is not crucial. 

In the inclusive cross section, the Standard Model physics background is larger 
than the signal. Furthermore, the detector-induced background from mismeasured jets in 
the forward region is several times larger than the real neutrino background, even for an 
ideal calorimeter covering r\ < 5.5.®^ First consider the lower-mass Cases I and II. Since 
gluinos and squarks are centrally produced with pT ~ M~, they give multiple jets and 
"round" events in addition to ^t- Jets with p^ > 75 GeV were found using the gemf ast 
fixed cone algorithm with R = 0.7. The minimum number of jets. A/jet, was varied between 
two and five. To identify round events, the sphericity in the transverse momentum plane, 
St, was calculated by summing all calorimeter cells with Et > 0.5 GeV and |?7| < 3. A 
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FIG. 47. signal for Case I MSSM parameters defined in Table 9 (open circles) and for 
QCD background (histogram) after requiring at least 5 jets with pt > 75 GeV and the 
sphericity and lepton veto cuts described in the text. 
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FIG. 48. Open circles: Ratio of signal and background curves from Fig. 47. Histogram: 
The same ratio for a perfect measurement. 

cut on St > 0.2 provided good separation of signal and background. After these cuts the 
signal to background ratio S/B for ~ 250 GeV was about 3 for Case I and about 5 
for Case II. The larger S/B for Case II reflects the harder spectrum from qji decays 
mentioned earlier. 



Semileptonic decays of gluinos and squarks are important; see Section 5.3.2 below. 
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FIG. 49. signal for Case II MSSM parameters defined in Table 9 (solid circles) and for 
QCD background (histogram) after requiring at least 2 jets with pT > 75 GeV and making 
the sphericity and lepton veto cuts described in the text. 

However, a lepton veto further improves the S/B for the distribution by rejecting ti 
and other Standard Model backgrounds. Events were vetoed if they contained a muon or 
an isolated electron. An electron was identified as an isolated electromagnetic cluster in 
the calorimeter with pT > 20 GeV and |?7| < 2.5, matched to a single track in the central 
tracker with a loose matching constraint, 

Av 

\E/p-l\ <max(0.5,3—) . (5.16) 

P 

Isolated muons with pT > 20 GeV and r] < 2.5 were identified using the standard gemf ast 
muon reconstruction. The efficiency of the lepton identification is not crucial for this 
analysis; even if it were perfect, there still would be background from r-decays of b and t 
quarks. 

The signal and background distributions for Case I with at least five jets and the 
sphericity and lepton veto cuts described above are shown in Fig. 47. The {S + B)/B ratio, 
shown in Fig. 48, reaches about 8 for = 250 GeV. Figure 48 also shows the {S + B) / B 
ratio obtained using flj^ calculated from the the missing u and Xi momenta, with the 
rest of the analysis unchanged. While the GEM calorimeter performance increases the 
background at low ^r, it is about as good as that of a perfect detector in the region for 
which the ratio is large. 

Figure 49 shows the signal and background for Case II, requiring at least two jets with 
PT > 75 GeV and the same sphericity and lepton veto cuts. For this case the direct decay 
Qr ~^ Xi9 dominates and leads to a significantly harder spectrum and to lower jet 
multiplicity. The {S + B) /B ratio is even larger in this case. Figure 50 plots the signals for 
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FIG. 50. Event numbers with > 250 GeV and the sphericity and lepton veto cuts 
described in the text vs. the minimum number A'jet of jets with > 75 GeV. Open 
circles: Case I signal. Solid circles: Case II signal. Histogram: QCD background. Case II 
has more events with low jet multiplicity because qr dominates. 

Cases I and II and the Standard Model background for > 250 GeV and St > 0.2 vs. 
the minimum number, A'jet, of jets with pT > 75 GeV. Both signals and backgrounds are 
constant for A'jet < 2 because it is impossible to have a large sphericity with only one jet. 
The signal falls off faster with increasing A'jet for Case II than for Case I because qr — > x'lq 
is dominant and gives a large rate for two jets plus ^t- Thus, the N^^t dependence provides 
a handle to distinguish among models. 

Given the large number of signal events, the statistical significance of the signals is not 
an issue. The t, W and Z backgrounds can be checked using isolated lepton samples; the b 
and c backgrounds can be checked using muons in jets. The resolution of the detector 
can be studied using inclusive data on QCD jets and on 7 + jets events. Given all these 
constraints, the background should be reliably known, so observation of a signal 5-10 times 
that expected from the Standard Model should be very convincing. The difficult problem 
of extracting the masses and other model parameters is briefly discussed in Section 5.3.3. 

For heavy gluino and squark masses such as those in Case III, is so large that the 
resolution is not important. Figure 51 shows the signal and background distributions 
for Case III with at least four jets having pT > 300 GeV and with the sphericity and 
lepton veto cuts identical to those for lighter masses. Heavy flavor backgrounds, detector- 
induced backgrounds from mismeasured QCD jets, and W and Z backgrounds are included. 
The QCD background dominates for low but falls more rapidly than the W and Z 
backgrounds, and so all must be considered. Since several hundred signal events survive 
these cuts with large S/B, it is evident that with an integrated luminosity of 10fb~^ GEM 
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FIG. 51. signal for Case III MSSM parameters defined in Table 9 after requiring at 
least 4 jets with pt > 300 GeV and making the sphericity and lepton veto cuts described 
in the text. Open circles: Signal. Solid squares: QCD background. Triangles: W ^ Iv 
and Z ^ vV backgrounds. Histogram: Sum of all backgrounds. 

could discover SUSY in this channel up to masses of order 2 TeV, about the upper limit if 
SUSY is related to electroweak symmetry breaking. 

5.5.2. Leptonic Signatures 

In addition to the IpT plus multijet signatures described above, there are many other 
signatures for super symmetry, including a number involving two or more leptons.®^ In 
particular, since the gluino is a self-conjugate Majorana fermion, 'gg and pairs can 
give isolated i^i^ pairs. Observing such like-sign pairs is essential for establishing the 
Majorana nature of any gluino signal. It also helps in separating gluinos and squarks. 
The dominant standard model background is expected to be from tt events in which 
either a 6 — > i.X lepton appears isolated or an isolated lepton sign is wrongly determined. 
These backgrounds, calculated previously in Ref. 86, are found to be negligible. For light 
gluinos, such as those in Cases I and II, the cross sections are so large that it is sufficient to 
use only the //^//^ channel, for which the lepton signs are very well determined in GEM. 
Therefore, only the issue of measuring the signs of electrons from Case III is addressed 
here. 

The same sample of Case III signal events described in the previous subsection was 
used for this analysis. While it is possible to enhance the leptonic sample by forcing a 
particular decay chain, e.g. — > XiQQ' ^ X\ ~^ Xi^^^: there are many such chains possible, 
no one of which obviously dominates. It was therefore decided to use the inclusive sample. 
For the background, only tt events, which are expected to dominate, were considered. A 
total sample of 30K tt events in ten bins with 50 < < 3200 GeV were generated, forcing 
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FIG. 52. pj^ distributions for the highest-py isolated lepton in dilepton events containing 
two isolated like-sign leptons. Open circles: Signal events generated with Case III MSSM 
parameters. Histogram: ti background from mismeasured electrons in the GEM central 
tracker. 

the decays t — > e^Ueb and i — > n^UiJ). This sample was used to determine the principal 
detector-induced background, that from misidentification of signs in the central tracker. 
From this, the total £^£^ background was determined. 

Electrons and muons with pT > 50 GeV and \r]\ < 2.5 were identified using the rela- 
tively loose cuts described in the previous subsection. These cuts, optimized for background 
rejection rather than for signal detection, appear to be adequate to identify this signal. At 
least two such leptons were required satisfying the isolation criterion 

Et <0.1pTe + 5GeV . (5.17) 

R=0.2 

Here, the prime on the sum indicates that the lepton itself is not included. This cut 
effectively rejects®^ the background from t i^vb and t bX^ b i^X. In addition 
a missing energy > 500 GeV and a transverse sphericity St > 0.2 were required. 
After these cuts, the total dilepton rates for the signal and for the ti background were 
comparable, so a very large rejection of unlike-sign pairs is not needed. 

Figure 52 shows the transverse momentum distribution of the highest pt lepton in like- 
sign dilepton events satisfying the above cuts. The lepton spectrum is soft relative to the 
gluino mass because the leptons arise from cascade decays. Figure 52 also shows the tt — > 
l'^vM~i>h background in which an electron sign is mismeasured by the GEM central tracker. 
The probability of mismeasurement was determined using the gemf ast parameterization of 
the central tracker electron resolution, including the tail from bremsstrahlung, as described 
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in Section 4.2. Muon signs are assumed perfectly determined, an excellent approximation 
at these momenta. Since the cuts not dependent on the electron sign reduce the background 
to the order of the signal, and since most of the signal Icptons have pT < 600 GeV, for 
which the central tracker determines signs with 95% reliability, it is not surprising that the 
background is small compared to the signal. The signal comprises several tens of events 
per 10fb~^ when both electrons and muons are combined. It would be uncomfortably 
small if one had to rely on only the n^/j.^ signal, which is a factor of four smaller. 

5.5.3. SUSY Parameter Determination 

In the MSSM there are at least eight mass parameters (M-, M~, Mj^, Mj , 

M~ , M~^, and Ma), two additional parameters related to the Higgs sector [fj, and tan/3), 
and still more parameters related to t decay. Nonminimal SUSY models have even more 
parameters. Since all supersymmetric particles in the MSSM ultimately decay into the x5) 
which is invisible, it is not possible to reconstruct any masses directly. 

The approximate mass scale can be inferred from the scale at which the signal 
deviates from the standard model background, as can be seen by comparing Figs. 47, 49, 
and 51. The mean ^or the distribution of the excess of events can be calculated very 
accurately for low masses because of the high statistics. However, the relationship of this 
mean to the masses is model dependent. For example, the missing energy is lower and the 
jet multiplicity is higher if M~> M~ than if M~ > M~. 

A large number of possible signatures may be used to determine the parameters. These 
include the cross section with multiple leptons,*^ multilepton cross sections arising from 
production of Xi^X2 ~^ similar channels, and the observed cross sections or limits for 

h,H 77; h,H ^ 4£] and t — ^ H^b. The tools to simulate these signatures have recently 
been developed,^^ and the methods to determine all of the MSSM parameters from these 
signatures are being studied. However, it is clear from the previous two subsections that a 
detector such as GEM, with reasonable calorimetry out to ~ 5.5 and tracking capabilities 
for charged leptons, is capable of observing clean samples of events in the relevant channels. 

5.6. High Mass Physics at Ultrahigh Luminosity 

While the studies presented above show that, at the SSC design energy of 40 TeV, the 
standard integrated luminosity of 10 fb~^ is sufficient for a wide range of interesting physics, 
there remain many potential new physics signals which require an integrated luminosity of 
order 100 fb^^. Practically speaking, this requires extended running at £ ~ 10"^^ cm^^ s^^. 
These new physics opportunities are almost all at very high mass scales, with effective 
cross sections of order 1 — lOfb. Examples include multi-TeV, weakly coupled Z' and W' 
bosons, quark-lepton substructure, and very massive technirho vector bosons. In these and 
most other examples, the new physics is signalled by an excess of well- isolated leptons at 
high pt. Robust and accurate measurement of such leptons was a primary design goal of 
GEM, as it is for the LHC detectors. 
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This section demonstrates the capacity of GEM at ultrahigh luminosity for precision 
measurements of the masses, widths and chiral couplings of 4 TeV Z' bosons via e^e~ and 
\i' final states and for the discovery and study of quark-lepton substructure at the scale 
A = 25 TeV in Drell-Yan production of high-mass dimuons. For both, we show that the 
GEM design permits distinguishing different models with a data sample of 100 fb^^. We 
also find that the reach of the GEM detector is M^/ ~ 8 TeV for new gauge bosons and 
A = 30 — 35 TeV for quark-lepton substructure. Finally, we briefly describe approaches 
to studying high-mass ^ — > i'^vi and quark-lepton substructure in qf^ — > ytr^v^x. We 
believe these are interesting because they effectively allow high-mass studies even when 
there is a neutrino and even when there is no well defined subsystem invariant mass.^^ 

The isolated high-energy lepton signals of these and most other new physics processes 
are relatively free of physics backgrounds. For example, a rejection factor of 7?.(e/jet) 
of C(10~^) reduces the jet backgrounds to Z' e+e~ to the level of a few percent. 
Straightforward isolation cuts on muons completely remove the main physics backgrounds 
to fi~^^~ . The most stringent requirements will be reliable operation of the trigger 

and detector and pattern recognition at such high luminosities. 

5.6.1. Precision Studies of New Heavy Z' Gauge Bosons 

Extensions of the standard electroweak gauge group, SU{2) involve additional 

neutral gauge bosons Z and sometimes charged W vector bosons. In the models of 
interest, the new weak bosons couple to quarks and leptons with strength of 0{e). In 
one year of running at £ = 10^^ cm~^ s"-*^, one could observe M^/ ~ 4 — 6 TeV, with 
10 Z' e^e~ events in a narrow Mg+g- range where none are expected. If such a 
boson were found, high-statistics studies would be needed to determine its nature. This 
subsection considers a 4 TeV Z' boson. For this mass, hundreds of events would be detected 
per year at ultrahigh luminosity. Since the couplings of such Z' bosons must be flavor- 
conserving, they may be detected in both Z' —>■ e+e" and Z' —>■ modes. Very 
precise measurements of the Z' mass and width can be made via the e+e" decay. As we 
shall describe below, the chiral nature of the Z' couplings to quarks and leptons can be 
investigated in the fJL^ jjT mode by measuring the distribution in the angle 6 between the 
outgoing iJL~ and the incoming quark. 

The goal of the studies presented here is to determine GEM's ability to distinguish 
two different Z' models by the bosons' widths and angular distributions. The two models 
considered were: 

1. A left-right model, in which SU{2) l ® SU{2)r ® U{1) breaks down to SU{2) l ® U{1). 
The SU {2)r coupling constant was taken to be the same as the SU (2) ^ one. The extra 
Z-boson of this model is called Zi below. 

2. A model in which the grand uniflcation group S'O(IO) breaks down to SU{5) (8) U{1), 
then to SU{2)l <8) U{1). Such a model may have an extra Z' = Z2 well below the 
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unification scale. 



The left-and right-handed couphngs to quarks and leptons in the two models are given by 



QuL 


= /9, 


9uR = 


/9 + 7; 


9dL 


= P, 


9dR = 




9uL = 


-3P, 


9vR = 


-3/3 + 7 


9lL = 


-3P, 


9lR = 


-3/3-7 



(5.18) 



The parameters /3 and 7 are 

(3 = 0.0528 , 7 = -0.3630 (Model 1) 
/3 = 0.0979 , 7 = -0.1958 (Model 2) . 



In calculating the Z' decay widths, it is assumed that there are three generations of 
quarks and leptons and that right-handed neutrinos exist and are much lighter than Mz'- 
If there are no other significant decay modes, the widths are given by 

nZ'^m = ^^^{9!L + s!K), (5.20) 

3 sm 20w 

where Cj = 3 for quarks and 1 for leptons. The full widths in the two models are r(Zi) = 
105 GeV and T{Z2) = 67.7 GeV. 

Z' e+e- 

PYTHIA 5.6 was used to generate 3000 events for each of the two Z' models, including 
full -y/Z/Z' interference.^^ Pileup at £ = 10^^ cm~^ s~^ was simulated by merging with the 
hard scattering process a set of minimum bias events. The latter were Poisson-distributed, 
with a mean number of 16 events per beam crossing. The total signal cross sections for 
3.75 TeV < M^+g- < 4.25 TeV are 5.44 fb for Model 1 and 4.54 fb for Model 2. No rapidity 
cut was imposed on the generated events. 

All tracks found by the central tracker (which, at 10^^cm~^s~^, has its inner silicon 
layers removed, leaving 8 planes of tracking) were processed and a most likely Z' vertex was 
found using momentum- weighted z locations. Then the 5x5 towers in the EM calorimeter 
were examined. Electron candidate towers were required to satisfy the following criteria: 

1. \r]\ < 2.46 so as to be within the tracker's coverage. Electrons falling in the transition 
region between barrel and endcap (1.01 < \r]\ < 1.16) were excluded because their 
energies are not as precisely measured. The geometrical acceptance of the detector for 
the Z' e+e- events is 78%. 

2. Et > 250 GeV in a 5 x 5 tower in the EC; 99.4% of the events with \r]\ < 2.46 passed 
this cut. 



78 



3. Transverse energy leakage into the 3x3 HC tower behind the 5x5 EC tower less than 
5% of the energy in the EC tower; 99.2% of the remaining events passed this selection. 
If a high energy EC cluster was isolated at this level, any energy in the first layer of 
the HC behind the cluster was added to the EC energy to improve the resolution. 

4. A charged track candidate with pT > 1 GeV must be within R = 0.0632 (Ary = 
0.06 X A(f) = 0.02) of the centroid of the EC tower. Essentially all electron candidates 
passed this cut. Much tighter track restrictions will be imposed below to reduce jet 
backgrounds. 

5. A calorimetric isolation cut of 



was applied. This accepted 99.6% of the signal events. 

The trigger efficiency for these events is close to 100%. Thus, GEM would collect 415 

Zi — i> e'^e~ and 350 Z2 e'^e~ events per 100 fb~^ run. 

Backgrounds to the Z' e^e~ signal come from misidentifying jets as electrons in 
QCD jet production and in ^ + jets with W ev. They may also come from 
isolated, high-mass e+e^ pairs from It production. The cross section for QCD jets with 
\(]\ < 2.46 and invariant mass Mjj — 4000 ± 250 GeV is 704 pb. The corresponding cross 
section for W + jets ^ e + jets is about 25 fb. The rate for tt — > e+e" + X with generated 
M-^^ > 3.5 TeV and pr(t J) > 100 GeV is 20.4 fb, for mt = 140 GeV. 

To be considered as background electron candidates, jets had to have Et > 200 GeV, 
one or two charged tracks in a cone of Arj x = 0.09 x 0.03 about the jet direction, and 
at least one charged track with pT > 50 GeV. This reduced a generated sample of 120K 
jets to 9360 (in 8890 events). These were then subjected to the cuts used to select electron 
candidates and 26 jets passed. This implies a jet rejection factor of 7?.(e/jet) ~ 2.9 x 10~^ 
at the 95% confidence level. The surviving dijet rate is 0.06 fb, less than 2% of the accepted 
signal rate. The VF+jet background is negligible. After the isolation cut, the tt e~^e~'+X 
background is completely eliminated by requiring Mg+g- > 3.5 TcV. (Sec the discussion 
below of the It fi^/j,^ +X background to Z' fj^^fjr.) When the requirement of exactly 
one reconstructed track with p^ > 50 GeV was applied to the Z' candidates, 83% of the 
events passed, for a net acceptance of 64%. 

The e^e~ mass resolution depends on the resolutions in both the electron energy and 
angle (see Eq. (5.3)). In the EC energy resolution, only the constant term matters for a 
4 TeV Z'. For the GEM baseline, this is 0.4%. (Any bremsstrahlung from the passage of the 
very high energy electrons through the CT material enters the same 5x5 EC tower which 
defines the electron shower and so does not affect the resolution.) For the high energies 
of interest here, the e"*" and e~ directions can be determined to an accuracy of about 




(5.21) 



i?=0.55 



5x5 
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2.4mrad using the longitudinal segmentation in the EC. This gives = 3.4mrad. The 
angular resolution can be significantly improved by using the tracker. The vertex location 
found using the two highest-i?r electrons has a resolution of 1.5 mm, while the EC hit 
location has a resolution of 0.25 mm for such high energy showers. The resulting angular 
resolution is = 1.5mrad. Then, since the typical opening angle between the electrons 
is greater than 60°, AMg+g- ~ 0.003 x M^/, dominated by the EC energy resolution. This 
is appreciably less than r^//2.35 for the two models under consideration. Indeed, models 
in which the the Z' width was artificially made as small as 25 GeV were considered and 
the width was shown to be resolvable.*^ 

Several luminosity-dependent effects were considered that might affect the Af^+g- res- 
olution.** Pileup noise in the calorimeters is unimportant. The effect of pileup tracks on 
the performance of the tracker for the purpose of vertex determination was included. 

Figures 53 and 54 show the Mg+g- distributions for the two Z' models as generated by 
PYTHIA and as reconstructed by gemf ast. The data in these figures correspond to about 
200 fb~^. These distributions are fit very well by a Lorentzian. The reconstructed masses 
determined from this fit are within 2 GeV of the generated 4TeV. The reconstructed 
widths are Vz^ — 123.0 ± 5.2 GeV and V = 88.5 ± 4.3 GeV. These are to be compared 
with the generated widths of 116.4 ± 5.0 GeV for Zi and 80.9 ± 3.9 GeV for ^2.*^ Within 
errors, the generated and reconstructed widths differ by an amount that corresponds to 
the detector resolution estimated above. This was verified for the artificially narrow Z' 
models. While the error in F^/ will be dominated by statistics, the main error in Mz' will 
be systematic, arising from possible nonlinearity in the EM energy scale in the TeV region. 

Finally, we estimated the reach of the GEM detector in M^/ for a 40 TeV collider 
with integrated luminosity 100 fb^^. Assuming that ten detected Z' — > e+e" events are 
sufficient for discovery (there is no background after cuts), the reach was found to be about 
8 TeV for both models considered. The reach for other models may be more or less by 
1—2 TeV. This high mass, then, was used to set the upper end of the dynamic range of 
the GEM electromagnetic calorimeter. 

Z' n+fi- 

The chiral couplings of Z' bosons to quarks and leptons determine the lepton angular 
distribution. Thus, the cos 9 distribution {9 is the angle between the incoming quark and 
outgoing fi~) can help distinguish alternative models. The angles and charges of high 
energy muons can be measured reliably even at ultrahigh luminosity. We imposed the 
criterion that muon charges must be known with at least 2.5a" significance. At least 99% 
of the signal events that passed the isolation and invariant mass cuts described below have 
the charges of both muons determined. Events with one mismeasured sign were rejected. 
The probability of two mismeasured signs is O(10~^) and no such events are expected in 
a data sample of 100 fb~^. 
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FIG. 53. The Mg+g- distribution for a 4TeV mass Z' of Model 1 as described in the 
text. Shown are the generated mass distribution (all events within acceptance) and the 
reconstructed distribution (passing geometric and calorimetric cuts 1-5). The number of 
events corresponds to 200 fb~^. 

To measure the cos^ distribution, one must determine, in addition to the muons' 
charges and directions, the direction of the incoming quark in the lab. The procedure 
usually followed in these studies has been to determine the (valence) quark direction by 
requiring that the Z' be produced at large rapidity, ?7^/ ^ 1.^^ Only about 20% of the 
heavy Z' —>■ jJi^ jJT events would pass this cut. We found it better to use all the events and 
plot the distribution in 

cos^* = sgn(r7B) tanh . (5.22) 

where rjE = ^{v^- + Vfi+) is the boost rapidity of the center-of-mass frame. ^'^'^^'^"'^ 

For > 2 TeV, the quark distribution is so much harder than the antiquark distribution 
that sgn(?75) determines the correct quark direction at least 75% of the time. Then 9* ^ 9, 
since the pT of the cm. system is generally much less than \/§. 

PYTHIA 5.6 was used to generate 1000 events of Z' with |M^+^- - Mz'\ < 

250 GeV for each of the two models described in the previous subsection. The cross sections 
in this simulation were found to be a{Zi — > ji^jjr) = 5.86 fb and a{Z2 — > jJi'^jJ^) = 4.49 fb. 
The geometrical acceptance of the muon system for all events was found to be 63% for 
both models; for muons with |?7| < 2.46, it is about 70%. For these muons, the trigger 
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FIG. 54. The Mg+g- distribution for a 4TeV mass Z' of Model 2 as described in the 
text. Shown are the generated mass distribution (all events within acceptance) and the 
reconstructed distribution (passing geometric and calorimetric cuts 1-5). The number of 
events corresponds to 200 fb~^. 

was highly efficient. Thus, before any additional losses, a data sample of 100 fb~^ yields 
375 Z\ — > n~ and 285 /^"""A*" events. A study of the muon system's geometrical 

acceptance for these Z' events was carried out using a full GEANT simulation, with results 
essentially identical to those obtained here with gemf ast. 

The only physics backgrounds to the Z' signal are continuum Drell-Yan production 
and ti [x^ \x~ + X. After momentum smearing is included, the signal region used for 
angular distribution studies was taken to be 3.2 TeV < M^+^- < 4.8 TeV. The Drell- 
Yan background was studied by generating 4000 'j/Z/Z' events with M^+^- > 1.5 TeV. 
Drell-Yan production contributed six background events to the signal region. 

The tt +X background was studied by generating 200K events with generated 

Mj^ > 1.4 TeV and ^^(t, t) > 100 GeV. The rate for this background is 1.4 pb, so the sample 
corresponds to 140 fb^^ (Other potential backgrounds, such as tW and WW production 
have very much smaller rates.) Detector response to events having M^+^- > 1 TeV and 
\rin\ < 2.6 at the generator level was simulated with gemf ast. There were 2900 such events 
(for 100 fb~^), of which 2700 passed the geometrical acceptance simulated in gemf ast. 
This is higher than the geometrical acceptance for Z' — > //''"/i" events because most of the 
muons from t-decay follow their parent's direction, which is fairly close to the beam, i.e.. 
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FIG. 55. The M^+^- distribution for the DreU-Yan continuum and Z2 H~^IJ- . The 
tt A*~'"At~ + X background before the invariant mass cut is also shown (dashed). 

away from the dead region near 77 = 0. Since these muons are usually accompanied by 
hadronic debris, the background can be removed by a combination of isolation and M^+^- 
cuts. 

An isolation cut using a large cone 



was made. Here, the sum is over energy in the full calorimeter and pt(a*) is the muon 
momentum corrected for energy loss in the calorimeter. This cut rejects 96.5% of the tt 
background while retaining 94% of the Z' events. All these events had identified muon 
charges. Finally, a cut on the measured dimuon mass of M^+^- > 3.2 TeV left no ti 
background events. 

A separate GEANT simulation showed that the rate of mismcasured low-py muons 
(e.g., ones with their tracks straightened by scattering) is small compared to the real high- 
ly muon rate. Therefore, the probability of getting two mismeasured high-py muons with 
Pt balanced to within 300 GeV is wholly negligible. 

A potentially serious loss to the Z' signal comes from muon-scattering debris 

that sprays into the muon system, degrading the reconstruction efficiency. This usually 
affects the first superlayer of the muon system. The average reconstruction efficiency for 
the high-energy Z' muons was found to be 85% per muon. the reconstruction efficiency 
for Z' — > events was raised to 96% by accepting events with at least one well- 

reconstructed muon plus one muon with sufficient hits in two unafi^ected superlayers to 
determine its angle. Since the pT of the poorly-reconstructed muon is expected to match 
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FIG. 56. Generated (dashed) and reconstructed (solid) cos 6* distributions for Zi — > 
events passing signal selections. 

that of the well-reconstructed one to within 300 GeV, its momentum can be determined 
with sufficient accuracy. If necessary to do so, it is reasonable to assume that the charge 
of the bad muon is opposite that of the good muon. 

In summary, the overall acceptance of the signal (including all produced events) is 
60%. There are 350 Zi and 270 Z2 n+n" events detected by GEM with 100 fb"^ of 
data. The Drell-Yan dimuon spectrum for the Z2 model is shown in Fig. 55, which shows 
all events that had generated invariant mass greater than 1.5 TeV and passed all cuts. 
For comparison, the distribution of dimuon events from the tt background that passed 
the isolation cut is also shown. The Drell-Yan spectra for the two models were fit with 
an exponential plus a Gaussian to determine the masses of the Z' enhancements. These 
were found to be Mz-^ = 4036 ± 40 GeV and = 3968 ± 65 GeV. The errors are purely 



The measured cos^* distributions for Models 1 and 2 are shown in Figs. 56 and 57. 
Events in the mass range 3.2—4.8 TeV were used for this analysis. There are 319 Zi and 257 
Z2 events in this sample; the missing events were lost due to momentum smearing. As noted 
above, only six Drell-Yan background events contaminate this sample. For comparison, 
the PYTHIA-generated distributions, corresponding to a perfect detector with GEM's ry- 
coverage, are shown in these figures. The agreement between theoretical expectation and 
simulated measurement is quite good. The forward-backward asymmetry is defined as 



statistical. 



Afb = 



N{cose* > 0) - N{cos9* < 0) 
7V(cos ^* > 0) + iV(cos ^* < 0) ■ 



(5.24) 
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FIG. 57. Generated (dashed) and Reconstructed (solid) cos^* distributions for Z2 
IJL^ IjT events passing signal selections. 

The results for the two models are 



The PYTHIA and gemf ast asymmetries agree well within the statistical errors and the 
central values of the two models are separated by about 3.5cr. 

5.6.2. Studies of Quark-Lepton Substructure in Drell-Yan Processes 

The greatest model-independent sensitivity to quark-lepton substructure is obtained by 
measuring the highest possible dilepton and dijet masses. If a signal is found, the inherent 
precision of lepton measurement at high energies will make the Drell-Yan process a much 
more incisive tool than high-£^y jet production for studying substructure. However, as with 
the Z' , hundreds of events at very large energies will be required. The contact interactions 
induced by quark-lepton substructure must be flavor-symmetric if the substructure scale 
A ^ 500 TeV. Thus, the Drell-Yan processes pp i~^£~ and i^Ui will have equal cross 
sections for electrons and muons. Signals in pp —>■ jJi^ jJT are considered here. While it is 
important to observe the substructure signal in pp e^e~ as well, the e^e~ spectrum 
cannot help distinguish between models because a difference in rate can be compensated 
for by a change in A. A brief discussion of substructure studies in the ji^v^, channel is 
presented in Section 5.6.3. 



Afb{Zi) 



0.113 ±0.056 (GEM) 
0.083 ±0.051 (PYTHIA) 

-0.090 ± 0.062 (GEM) 
-0.115 ± 0.058 (PYTHIA) 



Afb{Z2) 



(5.25) 
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Quark-lepton substructure modifies dimuon production in two ways: an excess of events 
at high M^+f^- = and a deviation in the angular distribution of the outgoing ji" 

relative to the incoming quark. ^"^ We consider two different chiral forms for the contact 
interaction arising from substructure at A = 25 TeV, well beyond that which can be 
studied at a 40 TeV hadron collider with 10'^'^ cm~^ s~^. The goal in this study is to 
distinguish these two models by their muon angular distribution. 

In the first model, only left-handed quarks and leptons are composite and the contact 
interaction is assumed to be the product of two weak-isoscalar currents (the "ISO" model): 



47r 



-^ISO = ~]y2 ^Lal'^QLa LLbl^LLb 



(5.26) 



where Qia = {iJ^a,da)L and L^^ — {va-i^aji are left-handed quark and lepton fields and 
a, 6 = 1,2,3 label generations. This interaction, together with the standard Drell-Yan 
processes, produces the subprocess cross section 



da{qiqi — > £ t^) 7ra 



2 r 



d{cos 9) 24s 
The functions Ai and Bi are given by 

4 



Ai{s) (1 + cosOf + Bi{s) (1 - cos 6*)^ 



(5.27) 



ah? 



sin^ 29w 



{T3i - Qi sin^ ew){i- sin^ ^^) 



s-M| 



+ 



Qi + Qi tan^ Ow 



Ml 



(5.28) 



Bi{s) 



+ 



Qi- 
Qi- 



1 



cos^ dw 

1 

cos^ 9w 



{Tsi - Qi sin^ 9w) 
Qi (i - sin2 9w) 



Ml 



s-Ml 



Here, sin^^^^ ^ 0.23, (Tg^, Q^) = (i, |) for qi = m and (-J, -J) for qi = di At high 
energies, s > ckA^, Ai{s) ~ {s/ah?)'^ and the angular distribution of l~ relative to qi is 
approximately (1 + cos^)^. 

In the second model the contact interaction is helicity nonconserving (the "HNC" 
model) : 

47r — 

-^HNC ~ ^ijQua'^Ra LLjb^Rb + h. C. , (5.29) 

where i,j 1,2 label indices in an electroweak doublet and eu = —621 = 1. This 
interaction, while theoretically unlikely, is studied here because it generates an angular 
distribution that becomes isotropic at large s. The interaction >Chnc affects UiUi — > £~£'^ 
only: 

da{uiUi — > £"£"*") TTO^ 



d{cjos 9) 



24s 



Au{s){l + cos 9Y + Bu{s){l - cos 9)' 
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+ 



TTS 

12A4 



(5.30) 
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Fig. 58. The M^+^- distributions for the ISO (a) and HNC (b) substructure models de- 
fined in the text. The lower dashed histogram is the standard Drell-Yan distribution. 

The contributions of >Chnc and do not interfere because their chiral structures are 

different. 

PYTHIA 5.5 was used to generate 1500 events each of standard DrcU-Yan (DY) and 
DY modified with the ISO and HNC contact interactions with A = 25 TeV. The muons 
were required to have generated Mu+^- > 2 TeV. The cross sections for events in which 
both muons have \-q\ < 2.46 are 2.83 fb, 9.44 fb and 7.88 fb for the DY, ISO and HNC 
cases. Detector response to events was simulated using gemf ast with Gaussian pileup for 
C = lO"^** cm^^ s~^. The acceptance was 70% for events with both muons in |?7| < 2.46, 
consistent with the acceptance of Z' — > events. Of these, 99.5% of the standard 

DY events had the charges of both muons determined. This dropped to 96% for the ISO 
and HNC models, refiecting the excess of high-energy muons generated by the contact 
interactions. 

The discussion of backgrounds and muon reconstruction is essentially the same as in the 
analysis of Z' The isolation criterion defined in Eq.(5.23) and the invariant mass 

cut M^+^- > 2.5 TeV eliminated the physics backgrounds. The acceptance of these two 
cuts for signal events passing previous selections was 41% for the DY events, 60% for the 
ISO events, and 66% for HNC. The higher acceptance in the ISO and HNC cases is due to 
the excess of high-mass dimuons. Retaining those events in which at least the momentum 
of one muon and the angle of the other are well-measured, the net reconstruction efficiency 
was found to be 95% for the DY and ISO cases and 92% for the HNC model. 
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Fig. 59. The cos 9* distributions for the ISO (a) and HNC (b) substructure models defined 
in the text. The lower dashed histogram is the standard Drell-Yan distribution. 

The net acceptances and number of events per 100 fb~^ are 23% and 80 events for DY, 
33% and 360 events for ISO, and 35% and 315 events for HNC. Figure 58 shows the M^+^- 
distributions for the ISO and HNC models detected in GEM per 100 fb"^ Also shown are 
the underlying DY mass distributions. The muon angular distributions were determined 
as in the Z' case. The cos 9* distribution for the fi~ in the ISO model, compared to DY, 
is shown in Fig. 59 (a) and for the HNC model in Fig. 59 (b). The DY background was 
not subtracted from the substructure-model distributions in these figures. The tendency 
for a (1 + cos^*)^ distribution in the ISO model and a flat one in the HNC model is clear 
and the two models are very well separated. The forward-backward asymmetries are 

AfbCDY) = 0.295 ±0.108 

Ai.B(ISO) = 0.328 ±0.050 (5.31) 
Ai7B(HNC) = 0.122 ± 0.056 . 
The errors are statistical only. Since the HNC and DrcU-Yan contributions do not interfere, 
the Drell-Yan contribution to the asymmetry can be subtracted. This gives ^^^^(HNC) = 
0.065 ± 0.065, consistent with an isotropic distribution. 

Finally, we estimate the reach in substructure scale A that could be attained with 
a data sample of 100 fb^^ at y/s — 40TeV. We define the reach to be that value of A 
which gives 5a over the Drell-Yan expectation, approximately 125 events. Then, since the 
subprocess cross sections above go as s/A^, it is possible to determine the reach by scaling 
the numbers of events found above. The reach of a 40 TeV pp collider estimated in this 
way is 30 - 35 TeV. 
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5.6.3. High Mass and Luminosity Physics Studies of Modes 

The range of physics that can be studied at ultrahigh luminosity can be greatly ex- 
tended by searching for isolated, high-pj- leptons accompanied by large — Pt- Ex- 
tensions of the standard gauge group generally involve W' as well as Z' bosons. Their 
mass and couplings can be determined in high-statistics studies of ^ — > e^i^e and iM^f^,-, 
respectively. With C = 10^^cm~^s~^ at a 40 TeV collider, it should be possible to carry 
out high-statistics studies up to M^/ ~ 5 TeV and to reach as high as 10 TeV. If contact 
interactions reflecting an underlying quark-lepton substructure exist, they may involve 
terms of the form udli' and its conjugate. Much can be learned about the chiral coupling 
of such interactions by comparing the rapidity distributions of the outgoing /j,'^ and 
High-statistics studies should again be possible up to A ~ 25 TeV. This section contains a 
brief description of the measurements that can be made at £ = lO^^ cm s if this new 
physics exists. 

Precise determination of the mass of a VF'-boson should be possible by measuring the 
Pt distribution of the electron in pp — > W' — > eu^. The main issue is how well the Jacobian 
peak determines M^r. For high M^/, there are no significant physics backgrounds. For 
example, the background jet -|- where the jet fakes an electron and Z*^ —>■ {—>■ Vu), is 
removed by the jet rejection 7^(e/jet) ~ 3 x 10^^. The detector-related issues are much 
the same as in the Z' e+e^ study, except that high precision is not as important here 
as the resolution on ^ 1 TeV. 

Information on the W' couplings to quarks and leptons can be obtained by measuring 
the angular (6) distribution in ^ — /U^z/^ decay. The ^-distribution of fi^ relative to the 
incoming d-quark will be the same as that of /x'^ relative to the incoming d. Thus, increased 
statistics can be obtained by adding data from both modes. These angular distributions 
provide the only way in pp colliders to detect the presence of right-handed neutrinos in 
the decay of W' bosons. To measure them, one proceeds as follows: (1) Select events 
with pt(a*) "near" the Jacobian peak found in W' — > eve and having balancing ^t- This 
enhances the signal relative to any Drell-Yan continuum or other background. (2) Measure 
the muon rapidity r]^. (3) Determine the neutrino rapidity, r]i,, by reconstructing the W'. 
For this, assume Pri^) = ^PTip)- As noted, this is a good approximation for multi-TeV 
W bosons since pt{W') < 200 GcV about 70% of the time for Mw ^ 2 TeV. The neutrino 
4-momentum and rjp are determined up to a quadratic ambiguity by imposing the W' mass 
constraint. One can either select the value that minimizes r]wri = 775 or accept each event 
twice. The analysis of the muon angular (cos 6*) distribution can now be carried exactly 
as was done for Z' and qq — > jJi^ jJT . Requiring good hits in all three muon superlayers, 
the fake high-py background to the W' — > -|- $t signal was found to be negligible. 

Contact interactions of the form udlv give an excess of muons at high-p^- Even though 
the parton cm. frame cannot be found in this case, it is still possible to obtain information 
on the chiral nature of the contact interaction by comparing the rapidity distributions, \'q^+ \ 
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Fig. 60. Iry^ I (solid) and \r]'J^ \ (dashed) distributions for the ISO (a) and HNC (b) com- 
positeness models defined in the text. 



and 1^7^-1 of muons with pT <^ 1.5 TeV. If, for example, the angular distribution between 
the incoming d-quark and the outgoing fj,^ is (1 + cos^)^ in the parton cm. frame, then 
|?7^-| is pushed to larger values because the d-quark is harder than the tt-quark and the 
fi~ tends to be produced forward. Correspondingly, the [r/^+l distribution is squeezed to 
smaller values. If the angular distribution is flat, as in the HNC model, the two rapidity- 
distributions will be identical. These features are illustrated for the ISO and HNC models 
in Figs. 60. These plots were made for the surviving dimuon event samples of Section 
5.6.2. The two models are well separated by the muons' rapidity distributions. If quarks 
and Icptons are composite structures, with a characteristic interaction scale in reach of 
hadron supercolliders, it is clear that measurements of the angular distributions in both 

contact interaction 



+ " and jj,^^ channels will be important for pinning down the chiral structure of the 
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6. CONCLUSIONS 

In this article we have presented simulations of the GEM detector's response for a 
selection of electroweak, flavor, and new gauge- inter action physics processes. The physics 
examples studied were chosen for their ability to help quantify GEM's capabilities for new 
physics and to aid in the detector's design and optimization, as well as for their intrinsic 
physics interest. While the examples form a fairly broad range of new physics within and 
beyond the standard model, they certainly do not represent a full survey of the interesting 
possibilities. 

The principal tool used for these studies was the GEANT-based parameterization of 
the GEM detector known as gemf ast. Although there were a few instances when there 
was no efficient substitute for full GEANT simulations of the response of limited portions 
of the detector — most notably in the calorimeter studies of jet backgrounds to H*^ — > 77 
— the gemf ast package was essential for our being able to carry out as broad a range of 
detailed simulations as we did in a fairly limited amount of time. We firmly believe it is 
necessary that the range of physics topics studied be as broad as possible. So long as the 
physics of the TeV-energy regime remains obscure, detectors built to explore it must have 
the capability to find whatever is going on within reach of their colliders. Thus, we believe 
that other large collider detector collaborations will profit from the development of similar 
tools for their large-scale simulation studies. 

Had the GEM project proceeded, further studies, both in depth and in breadth, would 
have been performed and an improved set of tools for modeling the detector developed. 
These in turn would have further influenced the evolution of the experiment. This iteration 
between design and physics simulation is an essential part of the creation of any modern 
detector. The simulations which were carried out guided the design of GEM and influenced 
a number of major technical choices. Examples are: 

• The requirements for the energy resolution of the electromagnetic calorimeter — both 
for the stochastic (~ 7%/\/E) and the constant (< 0.4%) terms — were set by the 
discovery potential for 77, H'^ ZZ* e+e^e+e", and Z' e+e". 

• The pointing capability and angular resolution of the electromagnetic calorimeter 
((40 — 50 mra d)/^E{GeV) + O.Smrad) were also motivated by the need to cleanly 
identify and measure H*-^ — > 77 and Z' — > e+e^. The pointing ability of the calo- 
rimeter complements the central tracker's determination of the event vertex, making 
robust GEM's precise invariant mass resolution for particles decaying into photons and 
electrons. 

• The segmentations of the calorimeter systems were dictated by the need to measure 
EM processes with high precision and the need to reject hadronic backgrounds. The 
calorimeter and the tracker combine to provide clean isolation for photons, electrons 
and muons. 
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• The precise momentum resolution of the muon system is needed to detect Z' — > pi^ii 
and H'^ ZZ* decays with high efficiency. 

• The resolution of the compact central tracker permits the complementary search for 
supersymmetry in the 'g'g — > + X channel. 

• The rapidity coverage of the forward calorimeters, extending to jr/l ~ 5.0, was dictated 
by the need to suppress backgrounds to the signatures for H.'^ — > t^i~Vv decays 
and for supersymmetry signals. 

• An adequate shielding system was designed to suppress neutron and photon back- 
grounds in the muon chambers, even at £ ~ 10^^cm~^s~^. 

While the simulation process was by no means complete, the results presented here demon- 
strate that GEM could have done outstanding physics at the SSC. We can only hope that 
others will have the opportunity to explore this physics elsewhere. 
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